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Quantum cosmology has traditionally been studied at the level of symmetry- 
reduced minisuperspace models, analyzing the behavior of wave functions. However, 
in the absence of a complete full setting of quantum gravity and detailed knowledge 
of specific properties of quantum states, it remained difficult to make testable pre- 
dictions. For quantum cosmology to be part of empirical science, it must allow for 
a systematic framework in which corrections to well-tested classical equations can 
be derived, with any ambiguities and ignorance sufficiently parameterized. As in 
particle and condensed-matter physics, a successful viewpoint is one of effective the- 
ories, adapted to specific issues one encounters in quantum cosmology. This review 
presents such an effective framework of quantum cosmology, taking into account, 
among other things, space-time structures, covariance, the problem of time and the 
anomaly issue 
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1 Introduction 



If quantum cosmology is ever to be part of empirical science, it must be described by a good 
effective theory There is no hope of exactly solving its equations in realistic models or to 
tame conceptual quantum issues made even more severe in the context of cosmology. The 
derivation of testable predictions requires systematic approximations at semiclassical order 
and beyond, and by experience with other areas of physics, effective actions or equations 
are the best available tools. 

While the speculative part of quantum cosmology, addressing for instance the Planck 
regime or the status of multiverses, requires all subtleties of quantum physics to be con- 
sidered — such as choices of Hilbert spaces, self-adjointness properties of Hamiltonians or 
unitarity of evolution, an understanding of deep conceptual issues and the measurement 
problem — physics can and must proceed without all these problems being solvedH Com- 
pare this situation for instance with quantum field theory, for which no rigorous interacting 
and non-integrable version is known. And yet, at the effective level it is the key tool be- 
hind the success of elementary particle physics. Given the immensity of the Planck scale, 
potentially observable effects in quantum cosmology are realized at low energies where 
semiclassical quantum gravity, with the first few orders in h taken into account, suffices. 
This feature makes effective theory in quantum gravity and cosmology even more powerful 
than in other settings [TJ [2]. As we will see in the course of this review, somewhat surpris- 
ingly, even conceptual problems of quantum gravity can advantageously be addressed with 
effective methods, especially with an extension to effective constraints. 

Given the amount of research on effective theories and their applications, one may think 
that deriving an effective theory of quantum cosmology is a simple and well-understood 
problem. However, this is not at all the case. Quite general and powerful techniques of 
effective actions or potentials do exist, employed with great success in particle physics 
and condensed-matter physics alike. Quantum cosmology, on the other hand, is unique by 
virtue of several features. It requires aspects of effective equations not encountered in other 
fields, related for instance to the prevalence of canonical methods, the generally covariant 
setting lacking evolution by a unique time parameter, or the absence or inapplicability 
of non-perturbative ground states or other distinguished classes of states. These technical 
problems will be discussed in due course. For now, as a motivation of our detailed look into 
effective theory, we state the following two general problems by which quantum cosmology 
differs from other fields. 

First, quantum cosmology is much like condensed- matter physics, with microscopic 
quantum degrees of freedom manifesting themselves on length scales far larger than their 
own. While the precise nature of microscopic degrees of freedom (strings, loops, . . . ) re- 
mains unclear and disputed among the different approaches, their presence in some form 

2 Especially in the loop-quantum-gravity community with its long and proud history of mathematical- 
physics primacy, it is sometimes said that one must become "less rigorous" in order to find useful and 
interesting physical results. This statement is, of course, incorrect; one does not become a physicist by 
being a mathematician first and then turning a little less rigorous. Physics requires as much rigor as 
mathematics, but a different kind of rigor. 
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is widely agreed upon. By considering large structures or space-time regions made from 
many building blocks, quantum cosmology differs significantly from most situations en- 
countered in elementary-particle physics, where events with a comparatively small number 
of particles in excited states close to the vacuum are studied. Quantum cosmology is a 
many-body problem, a situation in which it is difficult to derive and justify valid effective 
descriptions. The effective view has been put to good use in condensed-matter physics, 
but only thanks to rich and merciless experimental input to weed out wrong ideas and 
stimulate new successful ones. Quantum cosmology is not (yet?) subject to experimental 
pressure, and many (good and bad) ideas are sprouting. Effective cosmological theory 
must be able to stand on its own, requiring a systematic and rigorous formulation taking 
into account all features and consistency conditions to be imposed in quantum gravity. 

As the second problem, we observe that the theoretical foundation of quantum cosmol- 
ogy is much weaker than that of condensed-matter physics to which it is otherwise quite 
close. We know well which Hamiltonian we should use to find all states and energies of 
excitations in a crystal, but mathematically the problem is challenging and calls for the 
approximations of effective theory. In quantum cosmology, we don't even know which pre- 
cise Hamiltonian or other underlying object to use for the dynamics of a universe. Even 
if we choose one particular approach to quantum gravity, its mathematical objects or its 
specializations to cosmology are incompletely known or understood, opening wide the door 
for ambiguities and spurious constructions. 

We need a well-understood theory to pinpoint places where best to look for observa- 
tional effects, and we need observations to guide our theoretical constructions. Quantum 
cosmology, with an incompletely understood theory and no current observations, is a slip- 
pery subject, depriving us of a good handle to grasp its implications. In the absence of 
experiments, we can only rely on conceptual arguments and internal mathematical consis- 
tency conditions which, however, come along with their own problems. In this situation, 
effective techniques have proven to be one of the few reliable approaches for physical eval- 
uations of the theory, allowing one to include all crucial quantum effects in equations with 
clear physical meaning, and to take into account ambiguities by sufficiently general pa- 
rameterizations. This general framework and its current status in the context of quantum 
cosmology are the topics of this reviewjf] 

2 Cosmological consistency 

One source of technical and conceptual problems in quantum cosmology and quantum 
gravity in general is that the theory deals with relativistic space-time, in the absence of 
a unique Hamiltonian to generate evolution in time. Instead, many choices for time and 
corresponding Hamiltonians or evolution equations are possible, and they must all lead to 
the same physics. While this invariance is guaranteed classically, it implies a complicated 
problem after quantization, presenting the strongest set of consistency conditions to restrict 
possible choices of quantum cosmologies. Unfortunately, these issues are often set aside 

3 Fundamental issues of quantum cosmology are reviewed in the companion article [3]. 
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in research on quantum cosmology and even quantum gravity, owing to their complicated 
nature. This intentional oversight implies a large number of ambiguities, fixed in those 
contexts only by ad-hoc constructions. To keep this review focused, we will not discuss 
the rather large body of works in such directions, for instance those crucially using the 
distinction of a time variable by gauge-fixing or deparameterization in canonical settings, 
and only mention shortcomings in contexts in which they become apparent. Generally 
speaking, results derived with a distinguished choice of time (or gauge) cannot be consid- 
ered physically reliable unless one can make sure that they do not depend on one's choice 
of time. 

In this section, we will discuss the main features that an effective theory of quantum 
cosmology must deal with, which includes covariance and state properties. The former en- 
sures independence of choices of time, the latter deals with additional freedom in quantum 
theories. By these considerations, we will be guided toward suitable ingredients for the 
mathematical formulation of effective theory. 

2.1 Covariance 

Effective equations of quantum cosmology are supposed to modify some cosmological ver- 
sion of Einstein's equation by quantum corrections. In most cases, such a cosmological 
version is a reduction to isotropic Friedmann-Lemaitre-Robertson-Walker space-times or 
homogeneous Bianchi models (minisuperspaces), a restriction to some specific form of in- 
homogeneous degrees of freedom such as Lemaitre-Tolman-Bondi or Gowdy geometries 
(midisuperspaces), or an inclusion of unrestricted but perturbative inhomogeneity around 
a background in the former classes of models. 

2.1.1 Homogeneous models and automatic consistency 

In homogeneous models, the dynamics is completely determined by one equation for grav- 
itational degrees of freedom (such as the Friedmann equation) and one for matter (such as 
the continuity equation). The Friedmann equation of isotropic models, 



depends only on first-order derivatives and is therefore a constraint, to be satisfied by 
initial values of second-order equations of motion. When interpreted as a constraint (the 
Hamiltonian constraint), it is usually written in the form of an energy-balance law: 



with the matter energy -E ma tter — P^ 3 contained in some region of unit coordinate volume. 
(See [1] for a detailed discussion of coordinate factors when the volume is not fixed, espe- 
cially in the context of quantization.) In this form, the Hamiltonian constraint of gravity is 






(2) 
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obtained by varying the action by the lapse function iV = a/— goo- (For more on constraints 
and canonical gravity, see [5].) 

Given the Friedmann equation and the continuity equation 

p + 3-(p + P) = 0, (3) 

CL 

with pressure P, one can derive a second-order equation of motion by taking a time deriva- 
tive of ([1]) and eliminating p: the Raychaudhuri equation 

£ = + SF). (4) 

a 3 

At this stage, we have all equations expected from the components of the isotropic Ein- 
stein tensor: the time-time component providing the Friedmann equation, the (identical) 
diagonal components of the spatial part amounting to the Raychaudhuri equation, and all 
off-diagonal components vanishing identically. The equations obtained are automatically 
consistent with each other: By construction, the time derivative of the Friedmann equation 
vanishes if the Raychaudhuri equation holds. Therefore, if the constraint imposed by the 
Friedmann equation holds for initial values at some time, it holds at all times. 

This latter property is realized for a large class of systems describing versions of isotropic 
(or homogeneous) cosmology, not just for the classical one resulting from Einstein's equa- 
tion. As a phase-space function, H(a,p a ) in (T5]), with the momentum p a = —3(4irG)~ 1 aa as 
it follows from the variation <9Lg° av /<9d of the Einstein-Hilbert action reduced to isotropy, 
plays the role of the Hamiltonian generating all equations of motion in proper time. 
The Raychaudhuri equation indeed follows from the Hamiltonian equations of motion 
a = dH/dp a = {a,H} and p a = —dH/da = {p a ,H} (the first of which is identical to 
the definition of the momentum p a ). The matter terms p and P are realized by 

-^matter J E> ^-^matter /_\ 

P = ^T~ and F = ~^—fcT' (5) 

the negative change of energy by volume change. Matter equations of motion follow once 
-^matter in © is expressed in terms of canonical degrees of freedom. 

The phase-space function H(a,p a ) itself evolves according to the same Hamiltonian law, 
H(a,p a ) = {H(a,p a ), H} = and is automatically constant in time, no matter what form 
H(a,p a ) has. In homogeneous situations, the single Hamiltonian constraint that determines 
evolution is automatically preserved and consistent with evolution equations. We can easily 
modify H by any form of quantum corrections without encountering consistency problems, 
issuing a powerful license to cosmological model builders. 

Consistency remains valid when we consider different choices of time. So far, the 
equations were in proper time r. All other choices t in homogeneous models are related 
to r by r(t) = J* N(t')dt', with a lapse function N that enters Hamiltonian equations as 
well: If df(a,p a )/dr = {f,H}, we have 

df{a : Pa) = ^ d/( ?' Pa) = N{f, H} w { /, NH} 
dt dt dr J 
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for any other time with dr/dt = N. In the last step, we are allowed to pull the lapse 
function N inside the Poisson bracket, keeping the equality satisfied as a "weak" one, 
one that is valid provided the constraint H = holds. Applying the general law to the 
constraint itself, we again observe consistency: dH/dt = {H,NH} ps 0. 

The Hamiltonian constraint generates not only evolution with respect to a given time 
choice, but also the transition between different choices as a gauge transformation. In- 
finitesimally, with N close to one, we have r = t + e with e = f(N — l)dt, and for any 
function /, 

Sef = f(r) - f(t) =e^ = e{f, H} » {/, eH} . (6) 

The Friedmann equation, amounting to the Hamiltonian constraint, is automatically in- 
variant under changes of time, and so are its solutions. Also the evolution equations are 
invariant if we use the Jacobi identity for Poisson brackets: 

<l/ ,etf) = {{f,NH},eH} = {{f,eH},NH} + {f,{NH,eH}} 



dt 



d{f £ H} ~ {/, (de/dt)H} + {/, (5 e N)H} . 



Rearranging, we see that the gauge-transformed / evolves in agreement with the gauge- 
transformed df /dt, with a correction taking into account a possible gauge transformation 
of N and time dependence of e. 

The presence of a single Hamiltonian constraint therefore ensures dynamical consis- 
tency and invariance, even if quantum modifications occur. For this reason, homogeneous 
minisuperspace models are a simple and popular tool to investigate possible consequences 
of quantum gravity and cosmology. But for the very same reason, extreme care must be 
exercised when such models are used for physical predictions: The trivialization of consis- 
tency conditions does not hold in a more general context, and therefore spurious results 
can easily be produced in their absence. 



2.1.2 Inhomogeneity and covariance 

Compared with minisuperspace models, inhomogeneous cosmology presents a very different 
situation regarding consistency, even if inhomogeneity is small and treated perturbatively. 
For gravity and matter to fit together in Einstein's equation or a quantum modification 
thereof, a version of the contracted Bianchi identity must hold. But this identity, relat- 
ing different types of dynamical equations, is easily destroyed if quantum corrections, for 
instance those found in minisuperspace models, are inserted blindly. When amending in- 
homogeneous equations by quantum corrections, one must face the problem of anomaly 
freedom or covariance. Certain relations between different dynamical equations and gauge 
generators, classically implemented by the contracted Bianchi identity, must be preserved 
in the presence of quantum corrections. 

This well-known classical fact is often overlooked in quantum treatments, especially 
those making use of gauge fixing or deparameterization. Gauge generators disappear when 
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the gauge is fixed. A simple but illegitimate way out of difficult consistency problems is 
therefore to fix the gauge before quantum corrections are inserted. However, one then 
dispenses with ways to check consistency and cannot be sure that results obtained are 
physically viable. Most cases in which cosmological perturbations have been computed 
by gauge fixing in loop quantum cosmology, for instance, have by now been shown to be 
incorrect; important effects such as signature change have been overlooked. See also the 
instructive discussion of [6] in this context. We will come back to this issue later and for 
now continue with a classical discussion to provide more details. 

The first implication of the Bianchi identity is the presence of constraints. If we write 
V^G* = in the form 

<9oG° = -d a G^ - T V VK G* + r^G^ , (7) 

with spatial indices "a," it becomes evident that the components G° of the Einstein tensor 
cannot contain second-order time derivatives: On the right-hand side, all factors in the 
three terms are at most second order in time, and there is one explicit time derivative on 
the left-hand side, leaving only the option of first time derivatives in G° . These components 
of the Einstein tensor (minus 8ttG times the corresponding stress-energy components T° if 
there is matter) are constraints on initial values, while the remaining components provide 
evolution equations. In contrast to minisuperspace models, we are dealing with a larger 
constrained system of four independent and functional constraints, the Hamiltonian con- 
straint H = Gq — N a G° a and the diffeomorphism constraint D = NG®, which are to be 
imposed pointwise or for all possible multiplier functions N and N a in 

H[N) = - J d 3 xN(G° - N a G° a ) and D[N a ] = - J d 3 xN a NG° a . (8) 

We are dealing with an infinite number of constraints. (To define these integrations, one 
introduces a foliation of space-time into spatial surfaces t = const, also used to set up 
canonical variables. The time direction t a at each point, used to define time derivatives 
in evolution equations, may be different from the normal direction n a to spatial slices in 
space-time, a freedom parameterized as t a = Nn a + N a with the lapse function N and 
the shift vector field N a . The linear combinations of Gq and G° a in ([8]), depending on 
the shift N a and lapse N, take into account that constraints refer to directions normal 
and tangential to spatial slices, not to coordinate directions such as the zero-index of the 
Einstein tensor for a component along the time-evolution vector field t a .) 

As before, for consistency the constraints must always hold provided they are imposed 
for initial values, a feature that is guaranteed by the Bianchi identity as well. If we 
combine the Einstein tensor and the stress-energy tensor of matter in ([7]), we see that 
<9o(G° — 87rGT°) vanishes at any time provided the constraints G° — 8nGT® themselves (and 
therefore their spatial derivatives) vanish at that time and the evolution equations hold. 
By virtue of the contracted Bianchi identity, the constraints are consistent with evolution. 
Finally, as a third consequence, we see that all equations, Einstein's equation and the 
contracted Bianchi identity as a consistency condition, are covariant and independent of 
coordinates used. Therefore, solutions will be covariant. All equations hold irrespective 
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of the choice of coordinates, a well-known feature in perturbative cosmology which allows 
one to express all equations explicitly in terms of gauge-invariant variables [7J [8] . 

As in our discussion of minisuperspace models, a canonical view is useful to analyze 
which consistency conditions are satisfied automatically and which ones are non-trivial. 
There is a Hamiltonian constraint, and therefore Hamiltonian equations / = {/, C} 
are generated by a constraint C. However, in the inhomogeneous context, the con- 
straint is not unique (up to a pre-factor N), nor is time evolution. We have a much 
larger choice of possible time variables to generate evolution. The most general ver- 
sion of equations of motion is obtained if we use all our constraints in a linear combi- 
nation, defining H[N,N a ] := H[N] + D[N a }. For fixed N and N a , the Hamiltonian flow 
/ = {/, H[N, N a }} is then equivalent to Lie derivatives / = C t f along the time-evolution 
vector field t a = Nn a + N a in space-time, foliated by spatial slices with unit normals n a . 

For all constraints to be preserved by the evolution equations they generate, we need 
{H[M,M a },H[N,N a ]} ps for all M, N, M a and iV a . If this condition is satisfied, the 
constraints are said to form a first-class system. Unlike in homogeneous models, where 
H[M] always commutes weakly with H[N], the general condition is highly non-trivial and 
provides strong restrictions on consistent modifications of Einstein's equation. Quantum 
corrections can no longer be inserted at will. 

The same constraints that generate evolution provide gauge transformations, classically 
equivalent to coordinate changes. We use the same general combination as before, H[e, e a ], 
but interpret the multipliers e and e a differently, not related to a time-evolution vector field. 
Instead, the gauge transformation 8 e ^f — {/, H [e] + -D[e a ]} with the classical constraints is 
equivalent to a coordinate transformation or the Lie derivative C%f along the space-time 
vector field £ M with components such that e = iV^ and e a = £ a + iV a £° [9]. The factors of 
N and N a again result because space-time coordinate changes and the components of £^ 
refer to coordinate directions, while constraints refer to directions normal and tangential 
to spatial slices with normal n a = iV -1 (t a — iV a ). Also regarding gauge invariance, the 
condition of a first-class constraint algebra is then sufficient for consistency: In this case, 
all constraints are gauge invariant, 5 e H[N, N a ] = {H[N, N a ], H[e, e a ]} ~ 0, and so are the 
evolution equations they generate. With this full set of gauge transformations, one can 
freely change the constant-time spatial surfaces used to define canonical variables and to 
integrate the constraints ([5}. We are then dealing with a covariant theory of space-time, 
not just with a theory on a fixed spatial foliation. 

2.1.3 Hypersurface-deformation algebra 

The crucial consistency condition for any classical constrained system, its quantization, or 
an effective theory thereof, is therefore that it be first class: all constraints H[N, N a ] must 
have Poisson brackets that vanish when the constraints are imposed, 

{H[M, M a ] , H[N, N a ]} ps . (9) 
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For gravity, or any generally covariant space-time theory, the specific form is the hypersurface- 
deformation algebra [TO] 



{D[M a }, D[N a }} 
{H[M],D[N a }} 
{H[M],H[N}} 



D[C N bM a ] 
H[C N tM] 

D[q ab (MV b N - NV b M)] 



(10) 

(11) 
(12) 



with the spatial metric q . 

For a consistent quantization, an algebra of this form must be realized with commu- 
tators for constraint operators instead of Poisson brackets, and an effective constrained 
system must have quantum-corrected constraints such that a first-class algebra holds with 
Poisson brackets. If this is realized, no gauge transformations are broken by quantization 
and the quantum or effective theory is called anomaly-free. If this condition is satisfied, all 
consistency conditions that are classically implied by the contracted Bianchi identity hold 
in the presence of quantum corrections, and the (quantum) theory describes space-time 
rather than just a family of spatial slices. This property may be achieved with exactly 
the same form of the algebra, or with one that shows quantum corrections not just in the 
constraints but also in the structure functions of the algebra, as long as two constraints 
still commute up to another constraint. One universal example found in loop quantum 
gravity, as the most prominent result regarding hypersurface deformations with quantum 
corrections, has (fTUI) and fTTTT) unchanged, but (fT2|) modified to [UJ 



with some phase-space function f3. If the classical hypersurface-deformation algebra is 
modified, gauge transformations no longer correspond to Lie derivatives by space-time 
vector fields. Not just the dynamics but even the structure of space-time may be modified 
by quantum effects. We will discuss specific examples and results in later parts of this 
review. 

By analyzing quantum or effective constraints and their algebra, one can draw con- 
clusions about quantum space-time structures. For instance, once the full hypersurface- 
deformation algebra is known, one can specialize it to Poincare transformations by using 
linear iV and iV a . With N(x) = At + vx, for instance, we have a combination of a time 
translation by At and a boost by v. With /3 ^ 1 in f[T3"j) . the usual Poincare relations are 
modified. Although this may look like a version of deformed special relativity [T2 l [TB I [Tl]. 
there is no direct relation: In deformed special relativity, one has non-linear realizations 
of the Poincare algebra, with structure constants depending on the algebra generators. 
In (1T31 . we have corrections of structure functions depending on phase-space degrees of 
freedom, not directly on the algebra generators H[N] and D[N a }. A deformed version of 
special relativity would require a relation between phase-space variables, such as extrinsic 
curvature, and some space-time generators, such as energy. Relations of this form do exist 
in some regimes, for instance in asymptotically flat ones using the ADM energy, but not 
in general. 



{H[M}, H[N}} = D[f3q ab (MV b N - NV b M)} 



(13) 
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2.1.4 Consistent and inconsistent quantum cosmology 

In general terms, the problem of canonical quantum cosmology can be formulated as fol- 
lows: Find quantum corrections in H[N,N a ], perhaps motivated by some full theory of 
quantum gravity, such that these constraints remain first class. This statement presents 
a well-defined mathematical problem of classifying deformed algebras. For every consis- 
tent version that may exist, we can compute and analyze equations of motion by standard 
means, as explicitly written out above in the general classical case. At the present stage, 
several examples of consistent deformations of constraints remaining first class are known, 
mainly from loop quantum gravity, but there is no general classification of these infinite- 
dimensional algebras. 

As already mentioned, there are attempts to shortcut through the difficult calcula- 
tion of these algebraic structures by fixing the gauge before quantum corrections or other 
modifications are put in. If the gauge is fixed, one would no longer consider H[e,e a ] as 
gauge generators, but only use H[N, N a ] as constraints and to generate evolution. Some 
consistency conditions still need to be satisfied because the constraints are to be preserved 
by evolution, but this can usually be achieved more easily than in the non-gauge fixed 
case in which more fields are present. However, even if formally consistent versions of 
preserved constraints can be found in this way, such as those in [15J, they are not guar- 
anteed to be consistent because only a subclass of the constraint algebra can be tested 
when some modes are eliminated beforehand. Even if these are classical gauge modes, 
some consistency conditions and physical effects are overlooked. Moreover, the procedure 
is intrinsically inconsistent because one would first fix the gauge as it is determined by the 
classical constraints, and then proceed to modify the constraints that generate the gauge. 

A consistent scheme could be obtained only when the modified gauge structure is taken 
into account from the very beginning, but for that one would have to know a consistent 
version of non-gauge fixed constraints, not just of the gauge-fixed ones. Finally, having fixed 
the gauge, there is no way of calculating general gauge-invariant observables. Also here, one 
could only refer to the classical invariant variables, whose form however must be modified 
when quantum corrections are put into the constraints. Note that quantum space-time 
structures and modified constraints in most cases imply departures from classical manifold 
pictures, as we will see explicitly in the examples provided later. One can no longer refer to 
the usual form of coordinate transformations to compute gauge-invariant variables without 
using the constraints explicitly. In modified space-times, the constraints are the only means 
to compute gauge flows and invariants, but this can be done only if the gauge has not been 
fixed. 

Algebraic conditions can be simplified even more when gauge fixing is combined with de- 
parameterization. With the latter procedure, one chooses a phase-space degree of freedom, 
for instance from matter, to rewrite constraint equations as relational evolution equations 
with respect to this variable. Not just gauge transformations but even constraints then 
disappear from the system, and no strong consistency conditions remain. A popular ex- 
ample is the coupling of a free, massless scalar field 0, whose homogeneous mode in an 
expansion = + 5<p around some background can be treated as a global time function. 
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The canonical scalar Hamiltonian 



^scalar = \ J d*X f^== + V^^^) (14) 

+ f d 3 x8(---) J (15) 



1 ' - 2 



2 \ a/ det q 

with the momentum p^ = p^ + 5p^, then has a purely kinetic background term and is 
completely independent of 0. The momentum p^ is therefore a constant of motion and 
never becomes zero unless it vanishes identically. The background scalar has no turning 
points and is monotonic, serving as a global internal time along classical trajectories. 

The Hamiltonian generating evolution with respect to this variable is obtained by solv- 
ing the Hamiltonian constraint H griwity + H sca \ ar = (with fixed N as part of the gauge 

choice) for p^ = H^(q ab ,p ab \8<j),8p^): we have = {(f), H^} = 1 and p$ = {p^,,H^} = 0, 
consistent with evolution with respect to where the dot stands for d/d0, as well as Hamil- 
tonian equations for the remaining variables, such as q a b = {q a b,H^} and 8(f> = {8(f>,H^}. 
There is just one Hamiltonian generating all evolution equations, instead of a set of in- 
finitely many constraints. Almost as in minisuperspace models one can then implement 
in any quantum corrections one may desire. (In the context of cosmology, this ap- 
proach has been suggested for instance in [T6] in an analysis of Gowdy models and possible 
quantizations.) 

Deparameterization is a powerful mathematical tool to derive properties of physical 
Hilbert spaces, for which no general method exists in the absence of deparameterizability. 
(See also Section 14.1.51 ) In quantum cosmology, this method has been applied in [T7], and 
used recently to derive several aspects of self-adjointness and unitarity of evolution [18j [191 
BU] . But it cannot serve as a valid procedure to evade the anomaly problem and do physical 
evaluations. To start with, most realistic systems are not globally deparameterizable, with 
one variable having no turning points at all. Moreover, also this procedure, on its own, 
cannot weed out physical inconsistencies in spite of its formal consistency. One has the 
same drawbacks as in the gauge-fixed approach, and on top of that one has distinguished (or 
even introduced) one degree of freedom as time. For physical consistency, one should then 
show that results for observables do not depend on one's choice of time after quantization, 
but no systematic procedures exist to this end, constituting part of the problem of time in 
quantum gravity [2TJ [22j [23] . (See [21] for a discussion in the context of deparameterization 
or reduced phase-space quantization.) And even if one thinks that the distinguished time 
should be sufficient for all physical purposes and does not worry about relating results 
in different times, the fact that there is no analog of the Bianchi identity to constrain 
modifications should arouse suspicion. 

Gauge fixing and deparameterization are not necessarily bad, but they provide reliable 
results only when they are used after a consistent version of quantum-corrected constraints 
has been found. When this is the case, it is clear that all equations are consistent and gauge 
covariant, and instead of computing complete gauge-invariant variables for the corrected 
constraints, one may well pick a gauge or choose an internal time and work out physical 
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implications. There is also a possibility of providing consistent results with gauge fixing or 
deparameterization before quantization, but then one would have to show that all possible 
gauge fixings or deparameterizations would lead to the same physical observables. This 
requirement can be achieved in cases of simple constraints, such as the Gauss constraint of 
gravity in triad or connection variables as described at the end of Section 15.11 but presents 
a much more involved problem for the complicated Hamiltonian constraint. Problems and 
physical inconsistencies arise always when the gauge or time is chosen according to the 
classical system, and modifications of the constraints and thereby gauge transformations 
are inserted later. 

2.2 States 

Covariance or analogs of the contracted Bianchi identity severely constrain consistent ver- 
sions of quantum corrections, raising the manifold and metric structures underlying space- 
time in general relativity to the effective or quantum level. Although space-time structures 
and their covariance principles may then differ significantly from well-known classical ones, 
there is still a consistent dynamical theory independent of coordinates and gauge choices, 
and the set of all equations has meaningful solutions available for well-defined predic- 
tions. In addition to these restrictions on quantum corrections based on gauge aspects, 
the quantum theory of gravitational degrees of freedom itself should tell us what form of 
modifications we can have. 

2.2.1 Moments 

Effective equations, in general terms, describe quantum evolution by a smaller, more man- 
ageable number of degrees of freedom compared with the full quantum theory considered. 
When we go from classical physics to quantum physics, every degree of freedom we have 
at first is replaced by infinitely many parameters, for instance values the whole wave func- 
tion takes at all points in configuration space. It is difficult to deal with values of wave 
functions in physical terms. Another parameterization, more convenient for effective equa- 
tions as it turns out, is the set of expectation values of basic operators, (a) and (p a ) in 
Wheeler-DeWitt quantum cosmology, together with fluctuations and higher moments 

A(a b p c a ) := ((a - (a)) b (p a - (p a » c > Wey i (16) 

with operators in totally symmetric, or Weyl, ordering. For generic states, all these mo- 
ments with integer b and c such that b + c > 2, are independent of one another and of 
the expectation values, and therefore provide infinitely many degrees of freedom. (For 
b + c = 2, for instance, we have the two fluctuations A(a 2 ) = (Aa) 2 and A(p 2 a ) = (Ap a ) 2 , 
and the covariance A(ap a ) = \{ap a + p a a) — (a)(p a ).) 

Any set of effective equations provides a description between the classical limit and 
the full quantum theory of the system considered, making use of finitely many degrees of 
freedom for each classical one. There may be additional degrees of freedom compared with 
the classical theory, but not infinitely many more. By the new degrees of freedom included 
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and their coupling to expectation values — the variables with a direct classical analog - 
quantum corrections result. This is our general definition of effective theories, which we 
distinguish from coarse-grained theories. The latter may remove some quantum degrees of 
freedom by integrating them out, but typically leave whole towers of moments intact and 
therefore present infinitely many variables for some classical degrees of freedom. 

The behavior of the moments computed in simple unsqueezed Gaussian states with 
fluctuation parameter a [25] . 

A(a b p c a ) = 2 ~ ih+c)hC ° b ~ \ b/ 2yt/2)\ if b and c are even (17) 

and A(a b Pa) = otherwise, suggest a natural organization of effective theories of different 
degrees. A Gaussian state saturates the uncertainty relation 

{Aa) 2 {Ap a ) 2 -A{ap a ) 2 >^ (18) 

and, for position and momentum fluctuations of the same order, has a = Aa = 0(h}^ 2 ). 
The product h c a b ~ c ~ 0{h^ b+c ^ 2 ) then shows that moments of order n = b + c are of 
the order n/2 in h. We use this observation to generalize the notion of semiclassical 
states from simple Gaussians to a much larger class: If the moments of a state behave as 
A(a 6 Pa) ~ 0(h ( - b+c ^ 2 ) (said to obey a hierarchy), the state is called semiclassical. Moments 
of higher orders then affect only terms of high orders in ft and can be neglected in an 
approximation. Ignoring all moments of order higher than some n, only finitely many 
quantum degrees of freedom are left, and we obtain an effective theory. With n+1 moments 
of order n, an effective theory with moments up to order n has J2k=2(^ + -0 = | n2 + \ n ~ 2 
state parameters in addition to the two basic expectation values. 



2.2.2 Dynamics and quantum back-reaction 

Effective solutions, describing the evolution of a quantum state, depend not only on classical 
variables but also on the quantum state used, specified for instance by initial values of 
its moments. In quantum cosmology, one first promotes the constraint (J2]) expressed in 
canonical variables a and p a as 

H(a, Pa ) = — ^ + Matter 19 

6 a 

to an operator H, choosing some ordering of a and p a . Dirac quantization then implements 
the classical constraint equation H(a,p a ) = by the condition H\ip) = on physical states. 
In particular, the expectation value (H) must vanish in all physical states. If we express 
this equation as a functional equation on the space of expectation values and moments 
parameterizing states, we obtain an expression such as 
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where we have used (p 2 a ) = (pa) 2 + (Ap a ) 2 and the dots indicate additional terms that 
contain the covariance of a and p a as well as higher moments, and depend on the specific 
ordering chosen. If we take into account only expectation values in (H) = 0, the classi- 
cal Hamiltonian and the classical constraint surface are obtained. With fluctuations and 
higher moments, quantum corrections result that couple moments to expectation values 
and change the classical constraint surface. A systematic derivation and analysis of such 
moment terms, starting with expectation values of Hamiltonians and constraints, is the 
key ingredient of effective theories. 

The moments couple to expectation values, thereby providing quantum corrections to 
the classical motion, and their values therefore enter effective equations. A complete ef- 
fective theory cannot leave the moments in equations as unknowns and instead provides 
evolution equations or other conditions for them as well. But some freedom always remains, 
for instance in the initial values chosen for evolution equations of moments. If effective 
theories are formulated for expectation values without including additional degrees of free- 
dom such as moments to describe the evolution of classes of states, the specific states must 
either be restricted to provide unique effective equations, or be parameterized for a more 
general set. Often, such a state dependence enters effective theories only implicitly, for 
instance in the unique-looking low-energy effective action [26] free of any state parame- 
ters. This effective action describes low-energy effects of states near the vacuum of the 
interacting theory. In this way, the class of states is specified, certainly a rather small set 
compared to all possible states. 

In quantum cosmology, we are often interested in high-energy, Planckian phenomena 
and cannot restrict attention to the low-energy effective action. Even in low-energy regimes, 
which would be all we need to make contact with potential observations, it is not clear 
what low-energy state should be used. Quantum cosmology in its non-perturbative form, 
or quantum gravity in general, does not have a vacuum or other distinguished low-energy 
state. Effective theories of quantum cosmology must therefore be more general than the 
low-energy effective action, leaving more freedom for states parameterized in some suitable 
way. 

Even if we restrict attention to semiclassical regimes, the class of states to be considered 
may be large: simple Gaussians provide at most a 2-parameter family within a large set 
of semiclassical states when they are fully squeezed. For uncorrelated Gaussians, wave 
functions 



with a normalization constant N, we have, besides the two expectation values (a) and (p a ) 
on which the state is peaked, only one quantum parameter, the variance a. The most 
general Gaussian state is of the form ip z (a) = exp(— zia 2 + Z2d + z^) with three complex 
numbers z% such that Re^i > for normalizability. Out of these six real parameters, the two 
contained in z% do not matter for moments because the real part is fixed by normalization 
and the imaginary part contributes only a phase factor. For the remaining parameters, 




(21) 



15 



writing z\ = ol\ + and z 2 = a 2 + i/3 2 with real and fy, we compute expectation values 




{Pa) = ft- 



ctif3 2 - a 2 /3i 



(22) 



2«! 



and second-order moments 



(Aa) 2 



1 



(Ap a ) 2 = ft 2 ai + ft 




A(ap a ) = -ft 



0i 



(23) 



4«i 



2ai 



One easily confirms that the uncertainty relation ( !T8|) is always saturated. 

Deviations from saturation are much more difficult to parameterize, but easily occur 
for evolved semiclassical states even if they start out as a Gaussian. And even if one stays 
close to the saturation condition and does not vary second-order moments much beyond 
the values they can take for Gaussians, a Gaussian determines all higher moments in terms 
of the real and imaginary parts of z\. Our general semiclassicality condition A(a b p^) ~ 
0(ft( b+c ^ 2 ) provides an infinite-parameter family instead of the special 2-parameter one 
realized for Gaussians. Unless one can motivate Gaussians by other means, for instance 
proximity to the Gaussian harmonic-oscillator ground state or the vacuum of a free field 
theory, using them exclusively may easily be too restrictive. (Given the large parameter 
space of states, one may be tempted to refer to probabilistic arguments to pick "likely" 
states, following ideas that go back to an analysis of inflation in quantum cosmology [2T|,[28]. 
However, such probability considerations, though popular, are difficult, if not impossible, 
to make sense of in quantum cosmology [29].) 

Not just the absence of a vacuum state but several other special properties of quantum 
cosmology are important when we consider possible states: 

• Quantum cosmology considers long-term evolution. Even if we may be able to choose 
a specific form of semiclassical states at large volume and small curvature, it may 
change much when states are evolved back to high densities to infer possible impli- 
cations at the big bang. 

• As already mentioned, even the form of semiclassical states is unclear. It is customary 
to explore semiclassical features using simple and nicely peaked Gaussian states. 
In quantum mechanics, such states provide interesting information, and they are 
realized in exactly this form as coherent states or the ground state of the harmonic 
oscillator. In quantum field theory, Gaussian states are then close to the perturbative 
vacuum even for interacting theories. The dynamics of quantum cosmology, however, 
is not near that of the harmonic oscillator (except for some special models), and the 
unquestioned use of Gaussians is more difficult to justify. But going beyond Gaussians 
is complicated in terms of wave functions, whose parameters then become much less 
controlled. 

• It is not clear how precisely quantum cosmology can be derived from some full theory 
of quantum gravity, but the number of degrees of freedom is certainly reduced either 
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by exact symmetries or by perturbing around some background. In such situations, 
when degrees of freedom are eliminated, pure states easily become mixed. For general 
evolution equations able to model a full state, we should therefore allow for the 
possibility of mixed states, again going beyond pure Gaussians or other specific wave 
functions. Moments provide a parameterization of mixed states as well since they 
are based only on the notion of expectation values. 

• The question of covariance affects also the choice of classes of states. A state chosen 
in a minisuperspace model must have a chance of being the reduction of a full state 
that does not break covariance. This question may be difficult to analyze at the 
level of wave functions, but it also shows that a sufficiently large freedom in the 
choice of states must be included in considerations that aim to provide a consistent 
formulation of quantum cosmology. 

Since state properties are important for effective actions and quantum back-reaction, we 
should be as general as possible with the choice of states we consider. With wave functions 
or density matrices for mixed states, such a generality is difficult to achieve, but it is 
possible with parameterizations such as the one by moments. Moments provide a general 
form of semiclassical states, as already introduced. In effective theories, they are subject 
to their own evolution equations which show how they may change as high-density regimes 
are approached. They go well beyond the 2-parameter family of squeezed Gaussians, and 
describe pure and mixed states alike. 



2.2.3 Quantum phase space and covariance 

Covariance at the level of effective equations with quantum back-reaction can be addressed 
by combining the moment parameterization with the methods of the previous section, 
deriving consistent constrained systems amended by quantum corrections that include 
the moments. For the last question, we must be able to fit moments into a phase-space 
structure, so as to be able to compute Poisson brackets such as {H[M, M a ], H[N, N a }} 
with constraints that may include moment terms, such as an inhomogeneous version of 
(EH}. 

This construction is indeed possible: Together with expectation values, the moments 
form a quantum phase space with Poisson brackets defined by the commutator |30j, first 
for expectation values of arbitrary operators: 

{{A) t (B)} = ^p-. (24) 

This bracket satisfies the Jacobi identity and is linear. If we extend it to polynomials of 
expectation values by imposing the Leibniz rule, all laws for a Poisson bracket are satisfied, 
and we can apply the definition to moments. For instance, we obtain {(a), (p a )} = 1 and 
{(a), A(a b pl)} = = {(pa) , A(a b pl)} for all b and c. The moments, as defined here, 
are symplectically orthogonal to expectation values, a convenient feature for calculations. 
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Poisson brackets between different moments have been computed explicitly but are lengthy; 
see [30] and the correction of a typo in [25]. For low orders, it is usually more convenient 
to compute Poisson brackets directly from the definition (!24l) . For instance, we calculate 



{(Ag) 2 ,(A P ) 2 } 



{(q 2 )-(q) 2 ,(p 2 )-(P) 2 } 



(25) 




(HP 2 }) 



+ 4(g)(p) 



(I? 



ih 



p}) 



using the Leibniz identity and ( |24|) . 

If we know how a state enters effective constraints via its moments, a question which 
we will address in due course, we can compute Poisson brackets of effective constraints 
and see whether they provide a consistent deformation of the classical constraint algebra. 
If a consistent deformation is realized, the system can be analyzed further by standard 
canonical means to arrive at observables and dynamical equations. Note, however, that the 
Poisson tensor for moments truncated to some order is in general not invertible, for instance 
on the three second-order moments which form an odd-dimensional Poisson manifold: With 
( T2"5"|) and similar calculations for the other second-order moments, we obtain 



the 3-dimensional Poisson manifold of second-degree polynomials. Symplectic geometry 
therefore cannot be used for effective theories, while Poisson geometry is available from 
the definition (124)) . All relevant properties of constrained systems, such as the distinction 
between first and second class or properties of gauge transformations, can be formulated 
at the level of Poisson geometry [31] . 

2.2.4 Quantum-gravity states 

In summary of this section, we note that there are several specific issues in the derivation of 
effective theories for quantum cosmology, compared to other fields in which such methods 
are in use. The central theme is covariance, an issue which can be addressed only if one 
goes beyond the traditional quantum-cosmological realm of minisuperspace models. The 
notion of covariance and the question whether it is realized consistently may be affected by 
new quantum space-time structures introduced by the specific form of quantum geometry 
in one's approach to quantum gravity. 

The second, less obvious source of potential modifications of covariance is the form of 
dynamical quantum states used. Effective actions or equations depend on the classes of 
quantum states whose evolution they approximate, which should be specific solutions of 
some underlying quantum theory of gravity. Even if a theory such as quantum gravity is 
covariant, the selection of specific solutions may always break this symmetry. Covariance 
may then be realized in a deformed way, or only partially within one effective theory. 



{(Ag) 2 ,(Ap) 2 } 
{(Ag) 2 ,A(gp)} 
{(Ap) 2 ,A(qp)} 



4A(qp) , 
2(Ag) 2 , 
-2(Ap) 2 , 



(26) 
(27) 
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For instance, even if the states used are peaked on covariant classical fields, giving rise 
to covariant effective terms depending on expectation values, their fluctuations or higher 
moments may not be fully covariant. Changing the gauge in quantum gravity may then 
require the transition to a different effective theory, in which different state parameters have 
been used (reminiscent of the application of background-field methods in standard quantum 
field theory). On the other hand, if the class of states is sufficiently large, all gauge- 
related parameters could be encompassed within one effective theory. These considerations 
highlight the importance of the selection of states in the derivation of effective theories. 

Both quantum geometry and quantum dynamics must be part of one consistent quan- 
tum theory of gravity; in a complete treatment, their effects therefore cannot be separated 
from each other. However, they are derived by different means, making use of different 
expansions of the expectation-value function (H) of the Hamiltonian or Hamiltonian con- 
straint in the class of states used. In what follows, as in most derivations in the literature, 
we will split the treatment into one of quantum-geometry corrections first, as they are 
easier to see, followed by a discussion of quantum-dynamics corrections. When their ex- 
pressions are known, we can compare implications of different effects and see if some are 
more relevant than others in specific regimes. 

In the canonical setting, quantum-geometry corrections are specific to loop quantum 
gravity which has given rise to many results regarding background-independent quantiza- 
tion [32l E31 EI]. The effects in (H) are not unique, but characteristic enough to show 
implications for quantum space-time structure. After an overview of these terms in the 
next section, we will discuss dynamical quantum back-reaction, obtained from a further 
expansion of (H) by the moments parameterizing states. We will see the relation of canon- 
ical formulations to the low-energy effective action used in particle physics and the role 
of higher time derivatives in effective equations. Finally, we will put together our results 
to find properties of general effective actions taking into account all possible effects of 
canonical quantum gravity, and thereby shed light on quantum space-time structure. 

3 Quantum geometry of space 

Canonical gravity implements space-time structure by imposing the constraints H[N] + 
D[N a ] and requiring invariance under the gauge flow they generate. (There may be addi- 
tional constraints, such as the Gauss constraint if triad variables are used. Such constraints, 
however, restrict auxiliary degrees of freedom not related to transformations of space-time.) 
The diffeomorphism (and Gauss) constraint generates a simple flow by Lie derivatives, and 
has a direct action on quantum states. If we use states i^[q a b] as in Wheeler-DeWitt quan- 
tum gravity (leaving aside the complicated question of how to define an inner product) we 
have a formal action D[N a ]i/;[qbc\ = i>[£N a Qbc]- States are annihilated by the diffeomor- 
phism constraint if they depend only on spatial invariants, of the same form as classically, 
and expectation values in non-invariant states transform according to the classical gauge 
transformations. In loop quantum gravity, the action on states is different because discrete 
spatial structures do not allow infinitesimal diffeomorphisms to be represented. But for 
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finite diffeomorphisms, the situation is exactly as just described. Therefore, we do not 
expect significant quantum modifications of the spatial diffeomorphism constraint (but see 
Section l3~2l) . Similarly, the Gauss constraint, if there is one, has a simple and direct action 
that does not suggest modifications; see also the end of Section 15.11 

3.1 Hamiltonian constraint 

The Hamiltonian constraint H[N] presents a different story. It cannot be quantized directly 
by promoting its (complicated) classical gauge flow to an action on states. The only 
procedure to arrive at a constraint operator is tedious: inserting basic operators quantizing 
the classical kinematical phase space into the expression for the Hamiltonian constraints. 
As a phase-space function, the constraint contains rather involved combinations of the 
basic fields, which require regularization and sometimes even other modifications for a 
well-defined operator to result [351 EE]. Also on general grounds, we do expect quantum 
corrections in the Hamiltonian constraint because it includes all about the dynamics of the 
theory. The theory being interacting, quantum corrections must result. The Hamiltonian 
constraint is therefore the place where we should look for characteristic quantum corrections 
of different theories, as well as possible restrictions by consistency requirements. 

3.1.1 Gauge flow 

A constraint operator H[N] restricts states by H[N]\^jj) p h ys = and generates a gauge flow 
\ip) e = exp(—iH[e)/H)\ip). Physical states, on which exp(—iH[e)/h) acts as the identity, 
are gauge- invariant, and the gauge flow need not be considered separately if first-class 
constraint operators are solved. But exact solutions are complicated to find, and when 
quantum corrections are computed by systematic approximation schemes, the situation is 
quite different. Expectation values are constrained by (H[N]), a weaker condition than 
H[N]\ip) = 0. We have additional constraints ((/ — (f))H[N]) = 0, an enlarged set of 
constraints which all vanish in physical states, for arbitrary /. These constraints in gen- 
eral are all independent, constraining not just expectation values but also moments. In an 
effective theory, as spelled out in detail later, one solves this infinite set of quantum phase- 
space constraints order by order in the moments. To any given order, the sharp condition 
H[N]\i/j) phys = for physical states is not fully implemented, and on the corresponding 
solution space there is a non-trivial gauge flow by quantum constraints. Moreover, the 
effective treatment at this stage is more general because one does not assume a (kine- 
matical) Hilbert-space structure when solving the constraints for moments; therefore the 
standard argument of a trivial gauge flow does not apply. At the effective level, there are 
non-trivial gauge transformations by quantum constraints, bringing solution procedures 
closer to classical ones: one solves phase-space constraints and factors out their gauge. 

For an effective theory, the key ingredient is therefore the expectation value (H[N]) 
in a sufficiently large class of states, or the general expression parameterized by the mo- 
ments of states. As in (1201) . moments then appear in quantum corrections that change 
the constraint surface. With a modified constraint surface, the gauge flows must receive 
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quantum corrections as well for them to be tangential to the constraint surface mapping 
solutions to constraints into other solutions, as required for a first-class system. Indeed, 
the expectation value of the constraint, an expression that includes quantum corrections, 
determines the gauge flow on the quantum phase space with Poisson brackets (124j) by 

S e {6)(e) = ej- e (i>\eMiH[e)/h)deM-iH[e)/hM = ([ °'^ [c]]) = {(6), (H[e))} (29) 

for expectation values (0)(e) = t (ip\0\i(j) t , and therefore for any quantum phase-space 
function such as the moments by using the Leibniz rule. If the constraint surface changes 
by corrections in H, so do the gauge transformations. For the quantum corrected con- 
straint surface and the gauge flow we are therefore required to compute (H[e\) in general 
kinematical states. 

If quantum constraint operators are represented in an anomaly-free way, such that any 
commutator [Ci, is an operator of the form f(q,p)C3 with another constraint operator 
C 3 and structure functions f(q,p), the quantum constraints (Ci) form a consistent first- 
class system: {(Ci), (C2}} = ([Ci,C2])/ih = (f{.q,p)C 3 ). Since all expressions such as 
(f{q,p)C 3 ) vanish when computed in physical states, they provide effective constraints (in 
general independent of (C3) as phase-space functions). If all these constraints are imposed, 
the effective constrained system is first class and has a consistent gauge flow generated by 
all these constraints [23 EEl EH] • 

Such systems of infinitely many constraints for each local classical degree of freedom 
can be difficult to analyze. However, just as the moments obey a hierarchy in semiclassical 
regimes, the effective constraints can be truncated to finite sets to any given order in h or 
in the moments. The leading corrections can be found in direct expectation values (Ci), 
restricting expectation values with quantum corrections that depend on the moments. For 
instance, (1201) . when imposed as a constraint, shows how the classical constraint surface 
and the gauge flow of (a) and (p a ) are changed by fluctuations (Ap a ) 2 . These corrections 
depend on the value of (Ap a ) 2 , which is constrained and subject to gauge flows by higher- 
order constraints such as (p a H). When all constraints have been solved and gauge flows 
factored out to a certain order, the condition H\ip) = has been implemented for states 
used in expectation values and moments. One has then computed observables in physical 
states, sidestepping the complicated problem of computing an integral form of the physical 
inner product. This is one example for the use of effective methods to tame complicated 
technical and conceptual issues in quantum gravity. Even the problem of time can be 
solved, at least at the semiclassical level: Different choices of time are related by mere 
gauge transformations in the quantum phase space [101 SD 112] • We will come back to 
these conceptual question in Section 14.1.51 and for now note the lesson that expectation 
values of constraints supply the key ingredient for effective gauge theories. 

3.1.2 Quantum-geometry effects 

We should then apply effective techniques to the Hamiltonian constraint of gravity, com- 
puting (H[N]). Expressions for the Hamiltonian constraint are rather complicated even 
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classically, given by 



H^ab, * ) = --/= ^7T - -K) j + ^ (30) 

in ADM variables, the spatial metric q a b and its momentum 

lb7rG 

in terms of extrinsic curvature (with the spatial Ricci scalar ^R). 

Loop quantum gravity uses a different set of variables, the densitized triad Ef instead 
of the spatial metric q ab = E^E\j det(g c ^) and the Ashtekar-Barbero connection A\ = T l a + 
lK l a [131111], defined by combining the spin connection T l a compatible with the densitized 



triad and extrinsic curvature K % a = K a i,E bl / y/\ det(E^) |. The basic Poisson brackets are 

{Ai(x), E]{y)} = StttG^Cx, y) . (31) 

The Barbero-Immirzi parameter 7 > pHl H5] does not play a role classically but appears 
in quantum spectra and corrections. In these variables, the Hamiltonian constraint is 

tfgravK, E)) = - E J E }^— f (i4 + (1 + r 2 )e k m n(K - C)K - I?)) (32) 

with the curvature of the Ashtekar-Barbero connection. 

Both expressions are rather complicated to quantize, owing for instance to the ex- 
pressions of ^'R or T l a , to be written in terms of the metric or densitized-triad operators. 
Fortunately, there are several characteristic features, common to both versions: The Hamil- 
tonian constraint is quadratic in the connection or extrinsic curvature, and it requires an 
inverse of the spatial metric or the densitized triad. These two features, when combined 
with quantum-representation properties, imply characteristic structures of quantum geom- 
etry and associated corrections to classical equations. These corrections, in turn, can be 
analyzed for potential implications even if their precise form cannot be determined from a 
complete calculation of (H) in semiclassical states. And even if expectation values could 
be computed for a specific H and some states, ambiguities in the construction of H or 
the choice of states would be so severe that only general features and characteristic effects 
could be trusted. 

In a Wheeler-DeWitt quantization, we have formal quantum constraint equations with 
ir a b in ( 1301) replaced by functional derivatives with respect to q a b, acting on wave functions 
ip[q a b\- The formal nature leaves precise representation properties unclear, and therefore 
does not show specific quantum effects; one simply takes the classical expression and per- 
forms the usual substitution of momenta by derivative operators. One does not expect 
strong quantum-geometry effects, simply because quantum geometry has not been com- 
pleted in this setting. Loop quantum gravity, on the other hand, has celebrated its greatest 
success so far at this level of quantum representations, and indeed sheds considerable light 
on questions of quantum corrections resulting from quantum geometry. 
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3.1.3 Loop representation and background independence 

Loop quantum gravity looks closely at representation properties of basic operators. To 
eliminate the need for formal functional derivatives and the associated delta functions in 
quantizations of (EH), the classical fields A\ and are smeared or integrated over suitable 

sets in space: holonomies h e (A % a ) = V exp(J e d\t^A l a Ti) and fluxes Fef'(Ej) = J s d 2 yn^E^f l 
with the tangent vector t a e to a curve e and the co-normal nf to a surface S (on which 
an su(2)-valued smearing function f l is chosen). Allowing for all curves and surfaces, all 
information about the fields can be recovered. By the integrations, the delta function in 
{A l a (x), Ej(y)} = 87rjG6^6jS(x,y) is eliminated, and a well-defined holonomy-flux algebra 
results: 

{K{A% F S {E])} = 8 / K n fGh e -+ x Tif i (x ) h x ^. e (33) 

if there is only one intersection point {x} = eC\S, denoting by h e ^ x and h x -> e the holonomy 
along e up to x and starting at x, respectively. With strong uniqueness properties of 
possible quantum representations (1E1 El SHI EH] j spatial quantum geometry in this setting 
is under excellent control. 

Representations of the holonomy-flux algebra then provide operators to be inserted 
in f )32|) to obtain a quantized Hamiltonian constraint. The kinematical Hilbert space is 
spanned by cylindrical states ^f(A l a ) = ip(h ei ( y A t a ), . . . , h en (A l a )), each of which depends on 
the connection via a finite number of holonomies. The full state space, has no restriction 
on the number of holonomies that may appear, thereby representing the continuum theory 
rather than some lattice model. The inner product is obtained by integrating the product of 
two cylindrical functions over as many copies of SU(2) as there are non-trivial dependencies 
on holonomies in both states, using the normalized Haar measure. By completion of the 
space of cylindrical functions, the kinematical Hilbert space is obtained [50J. Holonomies 
h e then act as multiplication operators, changing the dependence of a cylindrical state on 
h e , or creating a new dependence if the state was independent of h e before acting. Flux 
operators, representing the densitized triad, become derivative operators in the connection 
representation used, and can be expressed in terms of invariant derivative operators on 
SU(2) (or angular-momentum operators). They have discrete spectra, indicating modifi- 
cations to the classical spatial structure j5B E21 E3] : distance, areas or volumes computed 
from El can, after quantization, no longer increase continuously even if the underlying 
curves, surfaces, or regions are deformed by homotopies. 

As a further consequence of discreteness, it turns out that holonomy operators do not 
continuously depend on the curves used. A cylindrical state depending only on h e and 
one depending only on h eoe ' with a piece e' appended to the curve refer to two different 
holonomies, and are orthogonal according to the inner product just described, even if e is 
a 1-point set. While we classically have t^,A l a {x)ri = ]im. e t^.M(h e >(A l a ) — l)/\e'\, where t a e , is 
the tangent vector of e' at x and \e'\ its coordinate length, the sequence h e i\ip) of states after 
quantization does not converge for any if e' is changed. The edge dependence disappears 
when one implements the diffeomorphism constraint, which has not been assumed in the 
previous constructions. States h e >\ip) for different e' are no longer orthogonal, but they all 
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give rise to the same state when diffeomorphism are factored out. The limit can then be 
taken, but always equals zero. Again, there is no way of deriving a connection operator. 

Unlike with classical expressions, it is not possible to take a derivative of h e (A l a ), for 
instance by the endpoint of e, to obtain an expression for a quantized A\. Loop quantum 
gravity does not offer connection operators; all connection dependence in the Hamilto- 
nian constraint must be expressed in terms of holonomies. No quadratic function as it 
appears in the constraint can exactly agree with a linear combination of exponentials, and 
modifications arise, motivated by background-independent quantum geometry. 

3.1.4 Holonomy corrections 

When the Hamiltonian constraint is quantized in loop quantum gravity, holonomies are 
used instead of connection components [35J, providing a "regularization" necessary to ren- 
der the constraint expressible by basic operators. However, the limit in which the "regu- 
lator" — the specific curves used for holonomies — is removed does not exist; we are not 
dealing with a proper regularization. (In [36], the limit is argued to exist and be trivial 
if spatially diffeomorphism-invariant states are used. But with this assumption there is 
no handle on the full off-shell algebra and its anomaly problem.) Instead, one usually in- 
terprets the difference of a holonomy- modified constraint or its effective versions with the 
classical expression as a series of higher-order corrections, amending the classical Hamil- 
tonian by higher powers of the connection, or intrinsic and extrinsic curvature. In this 
viewpoint, the modification is similar to expected higher-curvature corrections — except 
for the issue of covariance that we will have to address. (In some models of loop quantum 
cosmology, it is possible to represent an exponentiated version of the constraint in terms 
of holonomy operators without introducing modifications to the classical expression [54"] . 
However, the procedure seems to depend sensitively on specific properties of the model 
used and is not available in general.) 

Holonomy corrections therefore appear as higher-order terms such as 

- 2tr(rAK)) " = -^(t:K(xM4(x)MM*)) + Kt%A\{x) + ■■■ (34) 

obtained from a Taylor expansion of the exponential and a derivative expansion of the 
integration. These terms depend on the routing of the curve e, to be chosen suitably for a 
quantization of the Hamiltonian constraint. The condition of anomaly- freedom puts strong 
restrictions on the possible routings |55j which, however, are difficult to evaluate. Moreover, 
the expansion is done in an effective derivation, requiring the calculation of expectation 
values that lead to additional moment terms. All these calculations are difficult to perform 
explicitly, but the form of corrections is clear and quite characteristic: higher orders as well 
as higher spatial derivatives in the connection. With these types of corrections, suitably 
parameterized, one can look for consistent deformations of the constraint algebra to find 
versions for a physical evaluation of these corrections, or to derive further restrictions on 
the quantization choices made. 
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In effective equations, the use of holonomies as basic operators of the quantum theory 
has another consequence: The holonomy-fiux algebra (1331) is not canonical, with non- 
constant Poisson brackets. Moments of states used for effective descriptions are based on 
expectation values of basic operators, holonomies and fluxes in loop-quantized models. The 
general constructions of effective theories still go through: Poisson brackets on the quantum 
phase space, the quantum Hamiltonian or constraints, and so on. However, the explicit 
Poisson relations (1241) evaluated between individual moments of the form A(h b F c ) are 
different from the canonical case. In particular, moments no longer Poisson-commute with 
expectation values. This feature requires care and complicates some calculations. These 
problems, however, are only technical; see the examples provided later in Section 14.3.21 and 
in [561 [57]. 

3.1.5 Inverse-triad corrections 

Fluxes are linear in the densitized triad and do not suggest the same kind of modification 
as holonomies. Nevertheless, there is a characteristic effect associated also with them. Flux 
operators have discrete spectra, containing zero as an eigenvalue. Such operators do not 
have densely-defined inverses, and yet we need an inverse of the densitized triad for the 
classical Hamiltonian constraint. To obtain such quantizations, a more indirect route is 
taken in loop quantum gravity, which does result in well-defined operators but introduces 
another kind of correction, called inverse-triad correction. 

To see the form of these corrections, we show a lattice calculation of inverse-triad 
operators. As proposed in [361 [58] . the combination of triad components required for the 
Hamiltonian constraint (132]) is first rewritten as 

2^Ge^ k e abc ^£L= = Ul [ ^det^d 3 xl • (35) 
A/|det£/| 17 J 

In the new form, no inverse is needed, A l a can be expressed by holonomies, the volume 
operator can be used for J ^\ det E\d 3 x, and the Poisson bracket be turned into a com- 
mutator divided by ih. The resulting operators are rather contrived, especially with SU(2) 
holonomies and derivatives involved. But the presence of corrections and their qualita- 
tive form can be illustrated by a U(l)-calculation, for which we also assume regular cubic 
lattices. 

On a cubic lattice, we can assign a unique plaquette to each link e, by which we 
then label fluxes. Our basic operators are h e , which as a multiplication operator by 
exp(i f A a t^d\) takes values in U(l), and F e for a fixed set of edges e in a regular cubic 
lattice. We are therefore computing inverse-triad operators and their expectation values 
for a fixed subset of cylindrical states, but the lattice can be as fine as we want, allowing us 
to capture all continuum degrees of freedom. For this U(l)-simplification, the holonomy- 
fiux algebra, quantizing an Abelian version of (1331) . reads [h e , F e ] = —87fy£ph e while all 
operators commute if they belong to different edges. Moreover, the U(l)-valuedness of 
holonomies implies the reality condition h e h\ = 1, which we will make use of below. 



25 



We first rewrite ( 135]) in terms of holonomies instead of connection components, and 
express the volume V = J a/| det E\d 3 x by lattice fluxes ^/ 1 -^1-^2-^3! per vertex, with 
the three fluxes through plaquettes in all three directions around the vertex: t^{A a , V} = 
ih e {h~\ y/\F!F 2 F 3 \} or, more symmetrically, \i(h e {h~\ ^/\F x F 2 F 3 \}-h; 1 {h e , y/\FiF 2 F 3 \}), 
computes the Poisson bracket at the vertex. Since h e commutes with all but one of the Fj, 



we can focus on one of them, y |F e |. Keeping the power more general, the quantization of 
some inverse power of flux takes the form 

/,n,<^^ ht\F e \ r h e - h e \F e \ r ht t , . 

( F r_1 sgnF) e = el e 8 e ' el 6 =: h ■ 36 

167rGr 7 £p 

For any < r < 1 we quantize an inverse power of F but need not use any inverse in the 
commutator. 

Following (59] , we can now easily simplify these operators, if we observe the relations of 
the U(l)-holonomy-flux algebra, together with the reality condition. These relations imply 



hl\F e \ r h e =\F e + 8Tri%\ r , h e \F e \ r hl = \F e - 8tt 7 ^ 
such that 

j = |F e + 87r 7 4r-|F e -87r 7 4r (3?) 



Eigenvalues of this operator can easily be computed, with eigenstates equal to flux eigen- 
states [HOj EE]- All eigenvalues are finite, as required for a densely- defined operator, and 
show how the classical divergence of |F| r_1 at F = is cut off. For inverse-triad correc- 
tions in effective Hamiltonians, however, we need expectation values of I e in semiclassical 
states. Explicit calculations would require good knowledge of semiclassical wave functions 
or coherent states. 

For general effective equations it is sufficient, and even more useful, to perform a mo- 
ment expansion, keeping the specific state free and parameterized by moments. Staying at 
the expectation-value order of effective expressions, we have 

- = + ar^r - \(P) - g^r + moment teims . (38) 

167rGr 7 £p 

Already to this order we see characteristic corrections (depending on h via the Planck 
length). Inverse-triad corrections therefore have a contribution independent of quantum 
back-reaction. 

Interpreting (F) =: L 2 as the discrete quantum-gravity scale (the lattice spacing as 
measured by flux operators), we find the correction function 

ar(L) ' = w = (39) 

that will appear in an effective Hamiltonian constraint. To leading order in an expansion 
by h (or £p/L 2 ), the correction function equals one. But even if no moment terms are 
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included, there are quantum corrections in the full form of a(L). Corrections are strong 
for L 2 ~ Sir'yip or smaller, typically in the deep quantum regime, where a(L) drops to zero 
at L — 0. However, even for larger L, a r (L) is not identical to one and implies interesting 
corrections. 

In addition to corrections contained in (1391) and quantum back-reaction from moment 
terms, the flux dependence implies corrections from a derivative expansion of the integra- 
tions involved, as already seen for holonomies. Moreover, non-Abelian holonomies do not 
lead to exact cancellations in the substitution of h e {h~ l , V} for ^{A^, V} and rather imply 
additional higher-order corrections by powers of A l a [52] • As noted in the context of holon- 
omy corrections, such extra terms mix with higher-curvature corrections. The leading term 
in (I39j) . on the other hand, shows a different dependence on parameters that distinguish a 
given cosmological regime and are more characteristic. Their effects can thus be studied 
in isolation. 

3.2 Diffeomorphism constraint 

We have already stated that the diffeomorphism constraint can be quantized by its di- 
rect action on spatial functions or other objects such as curves and surfaces. In loop 
quantum gravity, for instance, a diffeomorphism $ acts by shifting all arguments of a 
cylindrical function by h e t— > /i$( e ), the usual pull-back of functions. The representation of 
the holonomy-flux algebra is diffeomorphism covariant under this action, showing that no 
quantum corrections to classical diffeomorphisms result. It is not possible to compute or 
represent an infinitesimal action or the diffeomorphism constraint because two states that 
differ by a non-trivial diffeomorphism are either identical (if the diffeomorphism does not 
change the underlying graph) or orthogonal. But finite diffeomorphisms suffice to remove 
the related gauge, which is done without quantum corrections. 

Nevertheless, the situation is not completely satisfactory because the diffeomorphism 
constraint is a crucial ingredient of the hypersurface-deformation algebra. If diffeomor- 
phisms are represented without quantum corrections, there should be no deformations 
of the relations (fTUj) and ( TTTj) of the hypersurface-deformation algebra for commutators 
involving at least one spatial deformation. However, the diffeomorphism constraint also 
appears on the right-hand side of (|T2|) . the crucial part for space-time structure. On the 
left-hand side, we have two Hamiltonian constraints, which we do quantize in loop quan- 
tum gravity and whose commutators we can, in principle, compute. The result should be 
a well-defined operator, which must vanish on physical states for the quantization to be 
anomaly-free. However, classically it corresponds to a diffeomorphism constraint, which 
cannot be represented directly. 

To check for anomaly freedom, one must then find an operator version of the right- 
hand side of ffl2|) . taking into account the structure function q ab , to be turned into an 
operator as well. This is one of the most important but still outstanding issues in loop 
quantum gravity, which was evaded by the arguments of and only partially addressed 
by the advanced constructions of [63j [M] . More recently, the issue has been revisited in 
several models [55] . with encouraging results. At least in U(l)-versions of 2 + 1-dimensional 
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gravity, one can indeed make sense of the right-hand side of (fl2j) as an operator, in such 
a way that the quantum constraint algebra is anomaly-free. As a side product, the same 
deformation ( Fl3l) with inverse-triad corrections as seen by effective methods [11] appears. 
(Holonomy corrections and their deformation of the constraint algebra could not be seen by 
the methods of [55], going back to [63j EI], because the consistency conditions of anomaly 
freedom are tested only at vertices.) 

In [55], the diffeomorphism constraint itself did not have to be amended by quantum 
corrections. However, other considerations in the same context have been put forward that 
may suggest such terms |65j. At present, the status regarding quantum corrections in the 
diffeomorphism constraint is incomplete, but a consistent implementation does not appear 
to be easy. From the point of view of effective theory, corrections to diffeomorphisms 
do not seem required because, in any canonical space-time theory, one is dealing with 
fields as functions on space. These functions are represented using some coordinates, but 
physics as always must be independent of the choice. There must therefore be a part of 
the gauge content of the theory that requires independence under arbitrary changes of 
spatial coordinates or, infinitesimally, invariance under spatial Lie derivatives. But then, a 
gauge transformation that amounts to a Lie derivative of all fields must have a generator 
identical to the diffeomorphism constraint uniquely associated with the fields [66]. The 
spatial structure assumed in canonical formulations leaves no room for corrections in the 
diffeomorphism constraint. 

The space-time structure is not presupposed in canonical quantum gravity and may well 
change, as indicated by some quantum corrections in the Hamiltonian constraint. Space- 
time, unlike space, has dynamical content and can easily receive quantum corrections, as 
borne out in loop quantum gravity. Having the classical structure of space but modified 
space-time is therefore consistent. Nevertheless, in an effort to relax some of the general 
assumptions of canonical formulations, one could expect changes to the spatial manifold 
structure as well, as perhaps indicated by potential corrections in the diffeomorphism 
constraint such as those in [65J. 

3.3 Quantum-geometry effects 

Comparing holonomy and inverse-triad corrections, we have several important properties: 

• Holonomy corrections crucially add higher powers of the connection to the classical 
quadratic form of the Hamiltonian constraint. In flat isotropic cosmological models, 
the connection is proportional to the Hubble parameter "H, which in turn is propor- 
tional to the square root of the energy density. Holonomy corrections in cosmological 
models therefore depend on the dimensionless parameters ipV, or a/ p/pp, both of 
which are tiny in observationally accessible regimes. 

Inverse-triad corrections, on the other hand, depend on the ratio £p / L 2 with the dis- 
crete quantum-gravity scale L in (13TJ1) . It is not easy to estimate L, but the dimension- 
less ratio associated with it certainly need not be small. Inverse-triad corrections can 
be more significant than holonomy corrections in observationally accessible regimes. 
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(The scale L may change in time, depending on the form of lattice refinement realized 
[63168].) 

• Holonomy corrections and inverse-triad corrections are both obtained from proper- 
ties of integrated objects, holonomies and fluxes. One should therefore expect not 
just higher-order terms as in the expansions already discussed, but also higher spa- 
tial derivatives in a derivative expansion of inhomogeneous models. For holonomy 
corrections, higher- derivative terms are crucial because they should be part of higher- 
curvature corrections together with higher powers of the connection that immediately 
arise from expanded holonomies. Only a suitable combination of higher powers and 
derivatives can result in consistent covariant versions. 

• Following up on the last item, we also need higher time derivatives to complete higher- 
order corrections to covariant objects related to curvature. Such corrections should 
be present even in homogeneous models, but are not easy to see directly from the 
form of holonomies. However, such terms cannot be ignored, because high-curvature 
regimes have significant contributions from higher-order and higher- derivative terms. 
In isotropic models, H 2 and Iri are of similar orders, both related to linear combi- 
nations of stress-energy components by the Friedmann and Raychaudhuri equations. 
An expansion of holonomies only by % (related to the isotropic connection; see Sec- 
tion I4.2p but ignoring higher time derivatives would be inconsistent. To see how 
higher time derivatives arise in canonical quantum theories, we have to pause our 
description of loop quantum gravity and return to more details of quantum back- 
reaction. 



4 Quantum back-reaction 

For a canonical effective theory, quantum Hamiltonians and quantum constraints (H), 
generating evolution or gauge flows by ( 129]) . must be expanded systematically by moments 
of states to see all quantum effects. This is also the case for individual non-linear correction 
functions such as (/) of inverse triads (138]) or (h) of holonomies as they may be implied 
by quantum-geometry effects of loop quantum gravity. Additional terms, products of 
expectation values and moments, are then added to the constraints. 

4.1 Effective quantum mechanics 

The correctness of the quantum dynamics resulting from a moment-expanded (H) can be 
illustrated with a quantum-mechanical example. We start with the well-known Ehrenfest 
equations 

ir = <»/•» • d -i = - {v ' m (40) 

for basic expectation values, computed using ( 1291 . These equations have been analyzed by 
[69] in the limit h — > to prove that quantum mechanics has the correct classical limit. 
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Effective equations go beyond this limit by performing a systematic expansion in H. 

The first Ehrenfest equation takes exactly the classical form, while the momentum 
expectation value is subject to quantum corrections: d(p)/dt = —(V'(q)) does not equal 
the classical force F((q)) = —V'((q)) at position (q) (unless the potential is at most 
quadratic). Moments as quantifiers of corrections arise when we expand the quantum force 
Fq = d(p)/dt as 

" (V'(q)) = -(V'((q) + (q- (q)))) = - £ ^ftfffW ) 

= F Q ((q),A(q n )), (41) 

or the quantum potential as 

V Q ((q), A(<f )) = (V((q) + (q- (<?)))> =V((q)) + jt ^^fW^ ( 42 ) 

n=2 U - 

such that -{V'{q)) = -dV Q /d(q). 
4.1.1 Quantum Hamiltonian 

The quantum potential is defined as a function on the infinite-dimensional quantum phase 
space of expectation values and moments, whose Poisson structure is given by (124]) . In a 
quantum Hamiltonian 

H Q = (H) = ^((P) 2 + A(p 2 )) + V Q ((q), A(q n )) , (43) 

we therefore have terms generating a dynamical flow of the moments by 

A(q b p c ) = {A(q b p c ),H Q }. (44) 

The coupled set of equations for expectation values and moments, (1401 and (1441) . is equiv- 
alent to the Schrodinger flow of quantum mechanics, but its solutions do not provide wave 
functions but rather variables directly related to observations. It can be solved with dif- 
ferent approximations, most importantly a semiclassical one by the order of moments, 
sometimes combined with an adiabatic one. In the latter case, applied to anharmonic os- 
cillators, effective equations are equivalent to those of the low-energy effective action [30] . 
The validity and usefulness of the canonical effective scheme is thereby established. 

In the context of quantum gravity, the feature of higher time derivatives in effective 
equations, a crucial ingredient of higher- curvature corrections, is of particular interest. The 
moments are related to such terms, although not in a direct way. Equation (HO"]) combined 
with fj4T]) already shows that a specific linear combination of the moments, with coefficients 
depending on expectation values, amounts to the time derivative of the momentum, or the 
second derivative of the position expectation value. Higher than second time derivatives 
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of (q) can be computed by taking further derivatives of (140]) and inserting A(q n ) and, 
for higher than third order, A(q b p c ) according to equations of motion (144]) generated by 
the quantum Hamiltonian. Different combinations of the moments therefore provide all 
higher time derivatives of (q). With the scheme just sketched, it is difficult to invert 
the equations to find expressions for moments in terms of higher time derivatives, or to 
eliminate all moments and end up with a higher-derivative equation just for (q) instead 
of the moment-coupled fHUj) . (J5TJ) and (PHI) . But with more-refined methods, as well as 
an adiabatic expansion, this task can be performed. For quantum cosmology, we learn 
that we must study quantum back-reaction to see all terms relevant for higher-curvature 
corrections. 

In semi classical regimes, the moments by definition obey the hierarchy A(q b p c ) ~ 
0(h^ b+c ^ 2 ), as can easily be verified in Gaussians; see (ITTj) and |25j. We can therefore 
consider the first term ^V"((q))(Aq) 2 for n = 1 in ( 142]) as the leading semiclassical correc- 
tion, providing a quantum Hamiltonian 

H Q = (H) = ^(P) 2 + V((q)) + ^(Ap) 2 + \v'\(q)){^) 2 • (45) 

(The kinetic term contributes (Ap) 2 /2m, potentially of the same order as ^V"((q))(Aq) 2 . 
But it does not appear in a product with expectation values and therefore does not cause 
quantum back-reaction.) For equations of motion of expectation values and second-order 
moments, relevant to this order, we use the Poisson brackets (126]) . Applied to our second- 
order quantum Hamiltonian, we find 



d(g) (P) 
dt m 



(46) 



(Apf-V"((q))(Aq) 2 (49) 
-2V"((q))A(qp). (50) 



(l/ V'((q})--V'"((q})(Aq) 2 (47) 

^ = 2 ~^P) (48) 
dt m 

dA(qp) 1 

dt m 
d(Ap) 2 

dt 

For a given potential, one may solve these equations numerically. However, it would be 
more instructive to compute (Aq) 2 and insert it in (f47|) to see what quantum corrections 
result. So far, all equations are coupled to one another and one cannot solve independently 
for (Aq) 2 (unless V" is constant, the case of the harmonic oscillator, a constant force or a 
free particle). But with an additional adiabatic approximation for the moments, decoupling 
can be achieved. 

4.1.2 Adiabatic approximation 

To zeroth order in an adiabatic approximation, we assume the moments (but not expec- 
tation values) to be time independent. We will denote the adiabatic order by an integer 
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subscript. Eqs. (HBl) and (!50|) then imply A (qp) = at zeroth adiabatic order, and (H9|) 
shows that (A p) 2 = mV A "((g))(A g) 2 . With the last equation, we see that the zeroth- 
order adiabatic approximation cannot be valid unless (q) is constant in time as well. To 
avoid such a restrictive condition, we proceed to higher adiabatic orders, from order i to 
order % + 1 by inserting time derivatives of Ai(q b p c ) on the left-hand sides of (HHl) -( l50l) to 
compute A i+ i(q b p c ) on the right-hand sides. (For a systematic implementation of the adi- 
abatic approximation, see [301 [70].) With time derivatives known from preceding orders, 
the equations to solve for the moments are initially algebraic, but additional consistency 
conditions relating different orders sometimes imply differential equations for coefficients, 
as we will see in this example. 
To first adiabatic order, 

AAqp) - V^) 2 - 1 d ( A °P) 2 (51) 

~ 2 m dt ~ 2V"((q)) dt [bL) 

1 /r(fd(g> 2 d[Aoqf\ 

using (HHl) . ( 1501) and our zeroth-order condition relating (Aop) 2 to (Aoq) 2 . The two lines 
can both hold only if 

d(Apg) 2 lV^M±d{q) 2 

= -2v^m^r {Aoq) ' (53) 

solved by 

(A g) 2 = ; ° (54) 

with a constant C. Our zeroth-order adiabatic relation between the moments then shows 
that (A p) 2 = mV"((q})(A q) 2 = mC^/V"((q}). Inserting these solutions in the quantum 
Hamiltonian f|45|) . we obtain a correction 

^-(A p) 2 + V((g»(A g) 2 = C^/WW) (55) 

to the classical Hamiltonian. 

As one goes to higher orders in the adiabatic approximation, one takes more and more 
time derivatives of A (q b p c ). We can see this feature already with the low-order equations 
found here. So far, we have used the first adiabatic order only to restrict the zeroth-order 
solutions. But with the solution (1541 found for (Ao^) 2 , we obtain from floTj) the moment 

AUav) - l m ^ q? - l Cm V "' m d(g) (56) 

depending on a first-order derivative of (q). We do not need A(qp) in the quantum Hamil- 
tonian, but this pattern continues for all moments at higher adiabatic orders, including 
Ai(q n ). When we go beyond second adiabatic order and insert solutions into expectation- 
value equations, higher-derivative effective equations will be obtained; see [70] for explicit 
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derivations. Quantum back-reaction by moments is responsible for these higher-derivative 
corrections, but there is no direct correspondence between the moments as independent 
quantum degrees of freedom and new degrees of freedom that appear in higher-derivative 
equations because more initial values need to be specified. It is not the moment expansion 
itself which gives rise to higher derivatives, but rather the adiabatic expansion of individ- 
ual moments. The order of moments corresponds to a semiclassical expansion, according 
to A(q b p c ) ~ 0(f^ b+c ^ 2 ) in semiclassical states, not to a derivative expansion. Any fixed 
order in h can produce arbitrarily high orders of time derivatives if the adiabatic expansion 
is pushed further. 

4.1.3 State dependence 

The parameter C in ( 154|) . related to second-order moments, is of the order h in semi- 
classical states; the correction ( )55|) is therefore the first-order semiclassical correction un- 
der the assumption of zeroth adiabatic order for the moments. We cannot choose ar- 
bitrary values for C because the uncertainty relation ( fl8l) must be obeyed, such that 
C = m _1,/2 AogAop > ^h/y/m. Requiring the uncertainty relation to be saturated deter- 
mines C. In general, this condition may be too strong because we would assume saturation 
at all times, amounting to the existence of a dynamical coherent state which is not guar- 
anteed for general potentials. But to zeroth adiabatic order, with the solutions found 
here, such an assumption is consistent: all dependence on (q) drops out in the product of 
(A g) 2 (A p) 2 (and we have A (qp) = 0). 

Without additional assumptions on the states solved for, or initial conditions for the 
moment equations (14"8]) - (|50]) . the constant C remains undetermined. One possibility to 
fix C, in the class of models of this example, is to assume that solutions are close to the 
harmonic-oscillator vacuum or some other specific state. If the potential V(q) = |mu 2 g 2 

is harmonic, (Ao^) 2 = C / \/ mu 2 is constant — in this case there are states for which the 
adiabatic approximation is exact — and equals the Gaussian spread a 2 in a coherent state: 
we may write C = o 2 \Jmoj 2 . For the harmonic oscillator, dynamical coherent states do 
exist and the uncertainty relation may be satisfied at all times. In this case, C = ^h/y/m, 
or (Ao^) 2 = ^h/mu, the correct relation for position fluctuations in the ground state. 
With (A p) 2 = \mftjjj, the non-classical terms in the quantum Hamiltonian amount to the 
zero-point energy ^hw. 

For a general potential V, we do not have the frequency parameter u to refer to, but 
we can define it as the square root of 2/m times the coefficient of the quadratic term 
in a Taylor expansion V(q) = Vo + V\q + \mu 2 q 2 + •••■, assuming that the coefficient 
is not zero. In this way, we treat higher than second-order terms in the potential as 
an anharmonicity. Specifying the class of states solved for as those that are close to a 
harmonic-oscillator ground state, we can therefore write (Aog) 2 = \J mV"((q}). The 
correction \h\fV'\(q)) /m in the effective Hamiltonian ( 1551) then agrees with that found 
for the low-energy effective action [26], a relation that holds to higher adiabatic orders as 
well [30]. 

The canonical picture of quantum back-reaction provides an interpretation of moment- 
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coupling terms as an analog of loop diagrams in quantum field theory, with moments 
taking the place of n-point functions. A formulation of the canonical effective scheme 
for quantum field theory is not fully worked out yet, but its implications for quantum 
gravity and cosmology can nevertheless be seen. Already in minisuperspace models there 
are characteristic effects which show cosmological implications of quantum corrections. 

4.1.4 Notes on the WKB approximation 

The WKB approximation is often seen as implementing a semiclassical regime, in the sense 
that leading terms in powers of h are considered in the quantum evolution equation for 
states, expanded as ip(q) = exp (ifir 1 Xl^Lo ^ n ^n{c[)) with an asymptotic series. With this 
ansatz in the Schrodinger equation, one can solve order by order in h to find expressions 
for the S n : in quantum mechanics, 



the classical Hamilton-Jacobi equation. 

Solutions obtained by the WKB approximation do not directly provide observables such 
as expectation values, for which additional integrations would be necessary. Such integra- 
tions are usually complicated to perform not just analytically but also numerically, given 
the strongly oscillating nature of WKB solutions in semiclassical regimes. Moreover, WKB 
solutions do not show how quantum corrections can be included in classical equations as 
the dominant quantum effects. In particular, although quantum back-reaction is implicitly 
contained in solutions to the WKB equations, it does not appear in the form of effective 
potentials or quantum forces useful for intuitive explanations of quantum effects. In the 
WKB approximation, So satisfies exactly the classical Hamilton-Jacobi equation, without 
any quantum corrections. Corrections to the dynamics arise by higher orders of S n in the 
wave function, but they do not appear in a form added to the Hamilton-Jacobi (or another 
classical) equation. 

While the WKB approximation, as an expansion in h, does have a semiclassical flavor, 
it can more generally be viewed as a formal expansion to produce solutions for wave 
functions. The WKB equations are obtained by solving the Schrodinger equation exactly 
at every order of h: An equation Yl^=o ^n^ n — is interpreted as implying E n = for all 
n. From a semiclassical perspective, on the other hand, one would interpret an equation 
J2^=o Enfr n = as providing a tower of quantum corrections Y^Li E n fo n to the classical 
expression E , and then be interested in solutions to the equations ^2^=0 E n fo n = cut off 
at finite orders of H. Additional consistency conditions are needed to determine the E n 
showing up in quantum corrections. Usually, the E n for n > depend on state parameters 
such as fluctuations, while E depends only on expectation values and equals the classical 
expression. A dynamical equation X^=o Enfo n = then encodes the quantum back-reaction 
of state parameters on the expectation values, implying deviations from classical behavior, 
as derived systematically by effective equations. 





for zeroth and first order in h implies S 
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In principle, one could derive such quantum corrections from WKB solutions by com- 
puting expectation values of the ^-expanded wave functions. But the WKB approximation 
does not automatically arrange the terms in its equation by semiclassical relevance. While 
canonical effective equations have a direct correspondence to the low-energy effective ac- 
tion, as already seen, the WKB approximation does not produce all terms [71]. Another 
question, important in the context of quantum gravity and quantum cosmology, is the 
treatment of quantum constraints (or the physical Hilbert space), which remains open 
in the context of WKB solutions. (For instance, one may solve H\ip) = with WKB 
techniques, but for approximate solutions, the gauge flow exp(— iH[e}/h)\ip)wKB does not 
automatically vanish.) Canonical effective techniques, on the other hand, apply to con- 
strained systems as well and even help to solve some long-standing conceptual problems of 
quantum gravity related to constraints and gauge. 

4.1.5 Effective constraints and the problem of time 

As already indicated in Section 13.1.11 a quantum constrained system with constraint oper- 
ators C produces quantum constraints Cq := (C), defined just like a quantum Hamiltonian 
( l4"3j) . but also independent quantum phase-space functions Cf = ((f(q,p) — (f(q,P)))C) (in 
this ordering) constrained to vanish in physical states [3T] 138] . In semiclassical expansions, 
calculating order by order in the moments, polynomial f(q,p) are sufficient. To fixed order 
in the moments, only finitely many constraints are then present. Their number is larger 
than the number of classical constraints because they remove not only expectation values 
of constrained degrees of freedom but also the corresponding moments. 

With the ordering of constraint operators to the right of f(q,p) chosen in effective con- 
straints, they are automatically first class if the constraint operators are first class. There 
are then constraint equations to be solved, and gauge flows to be factored out. The gauge 
flow is computed using the Poisson brackets (1241 . affecting also the moments. Standard 
techniques of constrained systems can then be used, except that moments truncated to a 
fixed order usually define a non-symplectic Poisson manifold. This feature requires some 
care and may have several consequences, for instance that the number of independent 
gauge flows does not equal the number of first-class constraints. Nevertheless, the usual 
classification of constraints and gauge flows applies |31j . 

To see the treatment of effective constraints we consider a Hamiltonian constraint op- 
erator C = p^—p 2 + W((j)) for a free, massless relativistic particle (q,p) coupled to a second 
degree of freedom ((p,p^) with an arbitrary ^-dependent potential W((p). Depending on 
the form of W(<f)), p^ may become zero along trajectories generated by the Hamiltonian 
constraint, in which case <fi does not serve as global internal time. On the other hand, 
with a g-independent Hamiltonian constraint, we could deparameterize by q, obtaining 

evolution by the classical Hamiltonian p = ±y^o| + W(<p). We have equations of motion 

d(p/dq = ±p 4> / yjp"^ + W{4>) and dp^/dq = t\W{4>)/ ^jp\ + W{4>). The momentum p^ 
evolves, and could indeed become zero. 
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The constraint operator gives rise to the effective constraints [101 |JT] 

C Q = (^) 2 -(p) 2 + (Ap ) 2 -(Ap) 2 + ^((0)) + iiy''((0))(A0) 2 (58) 

CV = 2(^)A(0^) + ^)-2pA(0p) + ^'((0))(A0) 2 (59) 

G n = 2(^)(A P0 ) 2 -2(p)A( w ) + iy'((0))(A(# )-^) (60) 

expanded to second order in the moments, together with additional constraints C q , C p , 
Cqp and so on, which we will not make use of. These constraints can be solved to find 
the quantum-corrected constraint surface, and their gauge flows can be computed to find 
moments of observables in physical states. Once the non-symplectic nature of the Poisson 
manifold of second-order moments is taken into account, these calculations are not very 
different from standard procedures. 

The effective constraints shown here illustrate another important feature: the complex- 
ity of constraints and their solutions. It comes about because effective constraints, to be 
first class, are defined in a non-symmetric ordering, while moments are by definition Weyl 
ordered. Reorderings required to express effective constraints as functions of the moments 
then introduce imaginary contributions by the commutator [q,p] = ih. For and C P4i 
to vanish, some moments must be complex. While moments before the imposition of con- 
straints, belonging to a kinematical Hilbert space, should be real as the expectation values 
of Weyl-ordered operators, after solving the constraints one moves to the physical Hilbert 
space, in general not related to a subspace of the kinematical one. After solving the con- 
straints, the original kinematical moments may therefore take complex values, as long as 
physical observables of the quantum constrained system are subject to reality conditions. 

For further consequences, we study the problem of time in this system, using the variable 
<p as internal time even though it does not deparameterize the system globally. At the full 
quantum level, local internal times, free of turning points only for finite ranges of evolution, 
cannot easily be made sense of: if internal time exists only for a finite range, evolution 
cannot be unitary even in this range. (See for instance the discussion in [721 E3 EI]. 
If states are evolved by local internal times past their turning points, evolution freezes: 
expectation values are stuck at constant values [751 EE]-) This consequence is the reason 
why the problem of time is much more severe at the quantum level, compared to the 
classical one. At the effective level, as we will see, the problem of time can be overcome, 
allowing consistent derivations of observables without using artificial deparameterizations 

gal mi p. 

If cf) is used as (local) internal time, it is not represented as an operator on the resulting 
physical Hilbert space, whatever it may be. No generally manageable techniques are known 
to derive physical Hilbert spaces and evolution in non-deparameterizable systems (for some 
possibilities of Hilbert-space derivations, see e.g. [751 ES E7J [78] ) . At the effective level, it 
is sufficient to distinguish as non-operator time by requiring that its moments in effective 
constraints vanish, 

(A0) 2 = A(0g) = A(#) = (61) 

while its expectation value {(f)) (denoted without the hat to indicate that <j) no longer acts as 
an operator) will become the time parameter. In fact, the conditions (IBTj) implement a good 
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gauge fixing of the second-order constraints C^, C p ,, C q and C p after quantization. (With 
constraints on a non-symplectic Poisson manifold, only three gauge-fixing conditions are 
required for four constraints. The remaining second-order moment involving <p, A((pp^), is 
fixed by the constraints, as we will see shortly.) In the terminology of [40J, these conditions 
implement the Zeitgeist during which ((f)) as local internal time is current. Imposing ( I6T1) 
initiates the transition to physical moments — moments computed for states in the physical 
Hilbert space on which does not act as an operator. 

Solving the effective constraints in the given Zeitgeist, we have A(0p^,) = —\ih from 

= 0, which then implies 

from C P4> = 0. Inserted in fl58l) . this implies the reduced constraint 

C = - (». + ^^(Apf + + m<0» (62) 

amounting to the quantum constraint Cq = (C) on the space on which and C P4> are 
solved in the given Zeitgeist. Solving C = for (p^), we obtain the time-dependent 
Hamiltonian for (0)-evolution, including quantum back-reaction. However, it still contains 
complex terms. 

In (|62|) . all terms except the last two are expected to be real- valued because (p) and 
Ap are physical observables, and (p^) can be interpreted physically as the local energy 
value. The constraint can then be satisfied, only if we allow for an imaginary part of (0), 
calculated from 

For semiclassical states, to which this approximation of effective constraints refers, we can 
Taylor expand the potential 

W((<f>)) = W(Re(<j>) + % lm(0)) = W(Re((f)}) + i lm(0) W'(Re((j))) + O((Im(0)) 2 ) 

by the imaginary term, expected to be at least of the order h because it vanishes classically. 
To this order, the imaginary contribution ImC = to C in ( 1621) implies that 

^) = ~y (64) 

The remaining terms, 

ReC = (p^) 2 - (p) 2 + M!_M!(Ap) 2 + H/(Re(0)) = (65) 

provide the physical Re (0) -Hamiltonian upon solving the constraint equation for (p^). At 
this stage, the Hamiltonian and its solutions, corresponding to evolving observables with 



37 



respect to 0, are all real: physical reality conditions are imposed and we have solutions 
corresponding to states in the physical Hilbert space. 

Although imaginary parts may be unexpected, a detailed analysis of this and other mod- 
els shows that they are fully consistent [H] . In models in which one can compute a physical 
Hilbert space, results equivalent with those shown here are obtained. Within the effective 
treatment of constraints, if we transform to a different internal time such as q, which is 
done by a gauge transformation in the effective constrained system so that a new Zeitgeist 
- the gauge-fixing f lBTj) — is realized, the imaginary parts are automatically transferred 
from ((f)) to {q), in such a way that observables remain real. By successive gauge trans- 
formations, one can evolve through turning points of local internal times, without freezing 
the evolution of physical observables; see in particular the cosmological example analyzed 
in [12]. The imaginary part of time can be seen as a remnant of non-unitarity problems of 
evolution in local-time quantum systems, but unlike in Hilbert-space treatments, it does 
not pose any problems at the effective level. 

Gauge transformations in effective constrained systems show that physical results are 
independent of the choice of (local) internal time. One may deparameterize the effective 
system in different ways to solve the resulting equations, without affecting observables. 
This conclusion, one example for effective solutions to the traditional problems of canoni- 
cal quantum gravity, indicates that deparameterization can be used consistently. However, 
in complicated systems subject to ambiguities such as factor-ordering choices, each depa- 
rameterization must be formulated in a specific way so that they all can result from one 
non-deparameterized system, effective or not. At the effective level, all quantum constraints 
and Zeitgeists must be computed and implemented with the same operator C for different 
local time choices to produce mutually consistent results. In many constructions of physi- 
cal Hilbert spaces, however, one quantizes a system with a specific deparameterization in 
mind, choosing factor orderings and using possible simplifications. In such a case, there 
is no guarantee that results can agree with those obtained from other parameterizations, 
and the independence of physical results of the choice of time is put at risk. 

4.2 Modified Friedmann equations, or: the sins of sines 

In quantum cosmology, a systematic derivation of quantum back-reaction is required es- 
pecially for reliable evaluations of holonomy corrections in the Hamiltonian constraint, as 
they both are relevant in high-curvature regimes and contribute to higher-curvature terms. 
Constraints appear in such systems, but for simplicity we will refer to deparameterized 
toy models. Holonomy corrections have provided a popular class of models within loop 
quantum cosmology, in which the Friedmann equation is modified in a simple way |79j : 
The classical constraint equivalent to the Friedmann equation is first modified to 



with a holonomy parameter i that could possibly depend on p. According to loop quantum 
cosmology [801 H] > t ne Hamiltonian is written in canonical triad and connection variables, 




(66) 
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with A\ = c5 % a and Ej = p5j, {c,p} = 8ttjG/3, under the assumption of isotropy. The 
densitized-triad component p can be positive and negative, according to the orientation of 
the triad, and is related to the scale factor by \p\ = a 2 . (Without loss of generality regarding 
effective equations, we take p to be positive in what follows.) For spatially flat models, 
as assumed here, c = 7a is proportional to the proper-time derivative of the scale factor. 
By the modification in fl66l) . the periodic form of holonomies is implemented, replacing the 
quadratic connection dependence of the classical expression. 

Computing Hamiltonian equations of motion for p allows us to eliminate c in favor of 
p, upon which the constraint equation takes the form of some kind of Friedmann equation. 
We have p = {p, H mo d} = (7^) _1 sm(2£c)y/p. Using trigonometric identities, we find 
sin 2 (£c) = \/l — 47 2 £ 2 d 2 ) with a = p/(2y/p). Inserting this in the modified constraint 
and solving for a 2 , the modified Friedmann equation becomes 

a) =— p { 1 -— ieap ) ■ (6?) 

This simple and interesting equation, with just a quadratic correction to the energy density, 
has served as the basis of many ad-hoc investigations of potential effects of loop quantum 
cosmology. 

Equation (167j) is clearly not an effective equation in the generality written here. Quan- 
tum back-reaction is ignored, while all terms in the complete series expansion of holonomy 
corrections 

sin2 ( fe ) _ 1 = + ±£V + . . . (68) 

{icf 3 + 45 + 1 ' 

are taken into account for the calculation. A consistent treatment would include only 
those higher-order terms in an expansion of holonomy modifications that are larger than 
any quantum back-reaction or other term that has been ignored. The relation of quantum 
back-reaction to higher-curvature corrections indicates that c 2 -corrections in ( 168]) should 
be of comparable size to c-corrections from quantum back-reaction, the latter of which are 
not included in ( )67|) . Including holonomy corrections but ignoring quantum back-reaction 
is therefore inconsistent, even at leading order in the c-expansion, unless one considers only 
models in which quantum back-reaction is weak. (Such models do indeed exist, as we will 
show later, but they are very special.) Keeping all terms in the c-expansion to arbitrary 
orders then leads to a questionable equation. 

One could think that keeping small higher-order terms is harmless, but it turns out that 
our cautionary considerations do matter for the form of modified Friedmann equations. To 
see this concretely, let us look at a few examples in which the holonomy modification 
in the Hamiltonian constraint is expanded first, followed by a calculation of Hamiltonian 
equations of motion and a modified Friedmann equation. The first order of c-corrections 
provides a constraint 
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Proceeding as before, we compute p = {p, Hi} = 27 1 c(l — ^£ 2 c 2 )^yp, solve for c in terms 
of p, insert the result in Hi, and rewrite as 



a\ 8nG ( 8nG 9 „ 9 9 . 

= —-p 1 - —iH 2 a 2 p . (70) 



Rather surprisingly, the result agrees exactly with the one obtained with the full holon- 
omy modification, ( 16 7p . However, this outcome does not mean that higher orders in the 
c-expansion do not matter. It rather shows that the specific form of the full modifica- 
tion by sin 2 (£c) is arranged so delicately that all higher-order contributions beyond the 
c 2 -correction precisely cancel one another. To confirm this, we go one order beyond the 
quadratic correction, modifiying the Hamiltonian constraint by 

H2 = - ( 1 - + + pp 3/2 = . (7i) 



8ttG7 2 V 3 45 

Again we proceed as before. (The higher higher-order polynomial equations to be solved 
for a relation of c to p can be handled easily within the perturbative scheme of the c- 
expansion.) The result, 

d\ 2 8-nG ( 8nG 9 „ 9 9 157 /87rG\ 2 a . a 9 \ 

-) = — p (1 - —yV* 2 P +«■(—) • (72) 

now has a higher-than-quadratic correction in the energy density. Going to higher orders 
in c and following this scheme shows that also the energy-order increases to include all 
possible powers. The leading corrections in p, such as (8nG/3)'y 2 £ 2 a 2 p with the same 
coefficient in all modified Friedmann equations, do not change if one goes to higher orders 
in the c-expansion and can therefore be used consistently — provided one stays in energy 
ranges in which it is the dominant term. When the energy density approaches Planckian 
levels and holonomy corrections are strong, however, 8nG^ 2 £ 2 a 2 p is close to one and all 
terms in the energy expansion are relevant. Bounce scenarios, for instance, cannot be 
formulated with a consistent version of the equation. (In this context, notice that the next 
term beyond p 2 enters with a positive sign. At high density, it may well be larger than the 
correction in (|67|) . in which case no zero of a and no bounce would be reached.) 

As anticipated, the sine-modification has its infinitely many higher-order terms ar- 
ranged such that all but the quadratic energy correction disappear. A consistent perturba- 
tive treatment keeping only the relevant orders instead produces a whole series expansion 
by the energy density. If one has reasons to trust the whole sine function and to exclude 
all other corrections, (|67|) is correct. But if there are additional corrections, however weak, 
it is not consistent to keep all terms in a c-expansion of the sine function; instead, one has 
a modified Friedmann equation with a perturbative expansion in p. Those additional cor- 
rections then unhinge the fine balance in the sine terms that eliminated all p-corrections 
beyond second order, and their form must be known for a reliable derivation of correct 
effective Friedmann equations. The main source of such extra terms is, of course, quantum 
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back-reaction, producing higher time derivatives that compete with higher powers of c. 
(Higher orders are also sensitive to quantization ambiguities, as analyzed for instance in 



4.3 Harmonic cosmology 

To understand the interrelation between different corrections, we should have a more de- 
tailed look at quantum back-reaction in quantum cosmology. In an effective description of 
Wheeler-DeWitt minisuperspace models, one considers the dynamics of expectation values 
(a) and (p a ) coupled to fluctuations and higher moments A(a b pl). The coupled dynam- 
ics, including quantum back-reaction, is usually complicated and unruly, but it simplifies 
considerably if perturbations around a simple model such as the harmonic oscillator in 
quantum mechanics can be used. As an analog of the harmonic oscillator in quantum me- 
chanics with simple effective equations, quantum cosmology has a harmonic model given 
by a free, massless scalar in a spatially flat isotropic geometry |56j . 

To realize the model as one with a Hamiltonian generating evolution, we must pick a 
time variable which we do by parameterizing, using the scalar <p as time. Since it is free 
and massless, the Hamiltonian constraint 



implies that p<§, is a constant of motion and has no turning points where p$ would 
move through zero. The scalar therefore provides a global internal time. Deparameterized 
models, as discussed before, cannot produce reliable physical predictions unless one can 
show that results do not depend on the choice of time. In the present context, we use 
the model merely to illustrate properties of quantum cosmological dynamics. For realistic 
effects, one can avoid deparameterization before quantization and the dependence on time 
choices by using effective constraints instead of effective deparameterized Hamiltonians 



It is an interesting coincidence that the same model is easily deparameterizable and at 
the same time, as we will see, harmonic, without quantum back-reaction. Both features 
imply that the model is extremely special even among symmetry-reduced isotropic systems; 
its implications must therefore be interpreted with a great amount of care. 

We perform deparameterization by solving the Hamiltonian constraint H(a,p a ,p^) = 

for 



the Hamiltonian generating evolution with respect to (p. Equations of motion for a(<f>) and 
p a {4>) are then obtained via Poisson brackets with p$(a,pa). If solutions are to be trans- 
ferred back to coordinate time, such as proper time, we solve d0/dr = {0, H^a^Pa.p^)} = 
p^/a(</>) 3 for 0(r) (with a constant p^) and insert this function in our solutions for a(</>) 
and p a {4>). 




(73) 



[33 EH na HQ. 




(74) 
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4.3.1 Effective Wheeler— DeWitt equations 

Effective deparameterized equations are generated by the quantum Hamiltonian (p^ (a, p a )), 

d{6) _ ([6,p^a,p a )}) 



d<fi ih 



{<0),<^(a,p a )>} (75) 



using the Poisson brackets f|24l . The absolute value inp^(a,p a ) makes a completely general 
expansion in moments complicated, but for [p^l not close to zero, there is a simple effective 
Hamiltonian. If we can ensure positivity of ap a in evolved states, the absolute value can be 
dropped. This is possible in particular for an initial state supported solely on the positive 
part of the spectrum of ap a (an operator for which we will assume Weyl ordering). Since 
ap a is preserved by the motion it generates, also after quantization, the evolved state will 
remain supported on the positive part of the spectrum of ap a . Unless |p^,| is close to zero, it 
is easy to find initial states supported only on the positive part of the spectrum of ap a and 
with specified initial expectation values for a and p a : Projecting out negative contributions 
will not change the basic expectation values much. We are then allowed to write (1751) as 



d + (0> + _ ± AixG + {[Q : \ap a \}) + ± l4nG + ([0,dpa})- 



d(f) V 3 ih V 3 ih 

= ± \fW {+{d) + ' +m + } (76) 

using |ap a ||'0) + = ap a \ip} + (and \ap a \ = dpj\ip) + ) on states with support only 
on the positive part of the spectrum of ap a . 

On such positively supported states, the 0-Hamiltonian is quadratic and can easily be 
expanded in moments. We have the quantum Hamiltonian 



H Q = ±^l^((a)(p a }+A(ap a )), (77) 

free of coupling terms of expectation values and moments: there is no quantum back- 
reaction. We compute and solve equations of motion for expectation values, resulting 
in 

(a)(0) = exp(±v/47rG/3 0) and (p a )(<f>) = exp( T ^4^/3 0) . (78) 



To transform to proper time, we solve d0/dr = p^ exp(=Fv / 127rG0) for 



= ± log(±Vl2^T) 



and obtain (a)(r) = (±V ll-nGp^r) 1 ^ , the classical dependence on proper time with a stiff 
matter source. (In particular, even after Wheeler-DeWitt quantization the system remains 
singular: infinite density p'i/(2(a) 3 ) is reached at finite proper time.) 
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In addition to these expectation- value solutions, the state evolves such that its second- 
order moments change by 



and dA(ap a )/d(p = 0, with solutions such that (Aa)/(d) and (Ap a )/(p a ) are constant. 
Semiclassicality is preserved exactly throughout evolution in this harmonic model, even at 
high density. Note that Aa and Ap a change nonetheless, but have often been assumed 
constant when state evolution was modeled by Gaussians. Wrong quantum corrections 
then result, which is especially significant in the presence of quantum back-reaction when 
the harmonic model is generalized. Especially curvature fluctuations are important because 
they grow when one evolves to high density, and they do show up in effective constraints 
such as ( 120]) as a simple example. This issue is another illustration of the importance of 
complete effective equations including the moment dynamics. 

With additional ingredients such as spatial curvature, a cosmological constant, a scalar 
mass or self-interaction, anisotropy or inhomogeneity, the system is no longer harmonic 
and becomes subject to quantum back-reaction. Deviations from the classical trajec- 
tory will then occur, to be captured by effective equations. (Some of these ingredi- 
ents also remove deparameterizability, but at the effective level we can still use local in- 
ternal times.) With a cosmological constant, for instance, the 0-Hamiltonian becomes 
P<t>( a ,Pa) = ±a/47tG/3 ay/pl — 4Aa 4 , a non-quadratic expression that entails coupling 
terms between expectation values and moments in a quantum Hamiltonian. 



4.3.2 Harmonic loop quantum cosmology 

At first sight, it seems that quantum-geometry corrections in loop quantum cosmology 
imply quantum back-reaction by deviations from the quadratic nature if fl66|) is used. This 
expectation is correct for inverse-triad corrections, but holonomy corrections, although 
they change the quadratic nature by higher-order terms, still lead to a harmonic model 
free of quantum back-reaction [56] . 

To see this, we first change to connection variables c = 7a = —(47r~fG/3)p a /a and 
p = a 2 , with {c,p} = 87T7G/3. The 0-Hamiltonian is still quadratic in these variables, 
proportional to \cp\, but the holonomy modification leads us to replace c by sm(£c)/£ with 
some £ that may depend on p. After this, the Hamiltonian p<j>(c,p) is no longer quadratic 
in c and p. However, if we introduce a new variable J := 3p exp(i£c) /87T7G , we have a 
linear 0-Hamiltonian 

p<t> = ±2 y -3 — — ( 81 ) 

if £ is constant. More generally, we can assume a power law for the p-dependence of 
£(p) = £op x , and define new basic variables 

V: = sJm-,) • (82) 
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so that the 0-Hamiltonian remains linear (and is just multiplied with 1 — x compared 
to ( )8T]) ). Moreover, and importantly, we have a (non-canonical) closed algebra of basic 
variables, 

{ V ,J} = -t£ J , {V,J} = i£ J , {J, J} = 2U V, (83) 

with the Hamiltonian a linear combination of the generators. 

Given these properties, upon quantization the Ehrenfest equations still provide closed 
equations for expectation values, without coupling to moments and quantum back-reaction. 
The Hamiltonian operator 

'j-jt 



±4/— -(l-x 



(84) 



is linear in J and its adjoint, and the quantum Hamiltonian is linear in ( j) and its complex 
conjugate. The only additional condition to impose is a reality condition because we have 
used partially complex variables. If we initially keep J and J as independent variables, 
valid solutions must satisfy J J = V 2 . 

We quantize by turning V and J into operators, choosing an ordering of J with the 
exponential to the right. The classical Poisson algebra is then replaced by the closed 
commutator algebra 

[V, J] = £ hJ , [V, J f ] = -l hP , [J, J f ] = -£ h(2V + £ h) , (85) 

where the last h comes from reordering exponentials. (Note that this non-canonical al- 
gebra implies (V, J)-moments not commuting with expectation values on the quantum 
phase space.) The reality condition, with the same ordering, reads JJ^ — V 2 = 0, which 
implies conditions on moments upon taking an expectation value, possibly preceded by 
multiplication with basic operators. For second-order moments, we have 

l(^)| 2 ~((V)+ t h/2) 2 = (AV) 2 - A(JJ) + \i 2 h 2 . (86) 

(For conditions on higher moments, see [82].) The reality condition is a Casimir of the 
commutator algebra of type sl(2, R) and therefore commutes with the quantum Hamilto- 
nian proportional to i( J — J'). If reality holds for initial expectation values and moments, 
it holds at all times. This statement is true not only for the harmonic model but for 
any Hamiltonian because the Casimir commutes with all V, J and Jt individually, and 
therefore with any function of these variables. 

Solutions for expectation values obtained from the linear quantum Hamiltonian are 

(V)((f)) = Aexp(C(f)) + Bexp(-C(f)) and (87) 
(J}(0) = Aexp(C0)-5exp(-C0) + ^p , (88) 

with two integration constants A and B as well as the constant C = ±2 A /47rG'/3(l — x). 
The reality condition then requires that 



\(J)\ 2 - (V) 2 = -4AB + ^ (89) 
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is of the order {V)£ h (the size of semiclassical fluctuations (AV) 2 and A(JJ) in (|86|) ). 
much smaller than p| for a universe which has a large amount of matter and is semiclassical 
at least once. (Recall that it is sufficient to impose the reality condition at just one time, 
for instance when semiclassicality is realized.) The product AB must therefore be positive, 
close to ilp^/C 2 , and the function (V)(<fi) oc cosh(C0— C<fio) never becomes zero. Holonomy 
modifications in the harmonic model replace the classical singularity by a bounce. 

As already seen in the beginning of this section, the bounce property of holonomy- 
modified equations is very sensitive to quantum back-reaction (or other possible correc- 
tions). With quantum back-reaction, the equations become more complicated but can still 
be analyzed numerically as long as moments do not become large. The approach to high 
densities can therefore be studied, but the Planck regime remains poorly controlled. In 
general, it is not known whether loop models always exhibit a bounce. 

4.3.3 Quantum Friedmann equation 

As in the case of holonomy-modified constraints, one can put effective equations into the 
form of quantum Friedmann equations. In the case of harmonic loop quantum cosmology, 

we eliminate (J) in favor of (V) by using the effective equations and, importantly, the 
reality condition. 

We have d(V r )/d0 = {(V), (p<f>)} oc Re(J) using f )84|) and fl85l) . which we turn into a 
proper-time derivative, as needed for a Friedmann-type equation, using d0/dr = p^,/a((j)) 3 
where a(0) is taken as some power of (V) (</>), depending on x according to (1821) . (More 
precisely, we write the Friedmann equation in terms of V by a/ a = V / (2V(l — x)) and then 
use (V)((fi). There is no direct relation between (a) and (V) because these are expectation 
values of different powers of the basic V unless x = 1/2.) The real part of (J), which is 
proportional to d(V)/d<j), is related to the imaginary part by the reality condition f )86|) . 
The imaginary part, finally, is proportional to the deparameterized quantum Hamiltonian 
p^. When all these relations are used and p^ is expressed via the energy density of the free 
massless scalar assumed here, we have 

Re(J) = ±(V r ) v /l-(87rG/3) 7 2 (^) 2 PQ (90) 

with 

Pq = Pfrcc + (3/87tG)( 7 ^)- 2 (A(JJ) - (AV) 2 )/(V) 2 (91) 

the free energy density corrected by moment terms. If (190!) is squared and suitable factors 
are inserted to express all terms by a, an equation for (a/ a) 2 follows: 

U =T fc ('-T*»4 (92) 

Except for the fluctuation terms in pg, this is Eq. ( 16T|) . 

There is another, more important difference: The derivation of (|9"2"|) holds only in the 
harmonic model, in which moments enter just by the reality condition, not by quantum 
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back-reaction. The quantum Friedmann equation is valid in this form only if the sole matter 
source is a free, massless scalar, and there is no spatial curvature, a cosmological constant, 
or deviations from isotropy. When any one of these conditions is violated, quantum back- 
reaction results and there are additional corrections, not just in pq but also in the general 
form of (1921) : As already anticipated by considering holonomy expansions in Section T4.21 
a whole series of corrections in a density expansion appears. The quantum Friedmann 
equation reads 

(a\ 2 8nG ( / 8ttG 2fe , 2 \ 1 F 8ttG ~ A2 w/ , a 6 W 2 2 \ 

UJ = -3- ^ - -3-^^^ ± 2 V 1 - -a-^"^^ ^ + "^-^ J (93) 

where W((f>) is a possible scalar potential, and we have a general quantum parameter 77 = 
r/fc + i(a 6 iy/p^) fc with coefficients rjk that depend on the moments, especially correlation 
parameters [831 El]- The expansion by a 6 iy//4 can be interpreted as one by (p—P)/ (p + P) 
with pressure P. (For a free, massless scalar, p = P.) A more-specific evaluation requires 
the detailed computation of quantum back-reaction to analyze how the moments evolve 
and what values they take especially at high density. Techniques for numerical studies 
have been provided in [82J. 

While it is difficult to find general information about the values of moments, it is clear 
that they contribute, among other effects, the canonical analog of higher-time derivatives 
as they appear in higher-curvature corrections. Moment terms and quantum back-reaction 
should therefore be large in high-density regimes, near the big bang. Reliable results 
in loop quantum cosmology can only state that the singularity is avoided by a bounce 
when matter is kinetic dominated, in which case W <C pkm and the ^-dependent terms in 
fl93|) can be ignored unless 77 is extremely large. This conclusion coincides with numerical 
investigations [851 ESI EZ] of the underlying difference equation for wave functions. However, 
such numerical studies suffer from the choices required for wave functions, for instance by an 
initial state. With such methods, it is difficult to capture general-enough effects, which in 
quantum cosmology with its lack of distinguished states make robust conclusions difficult. 
As shown by the generality of the quantum Friedmann equations displayed here, effective 
techniques allow one to draw conclusions and confirm the regime-dependent validity of 
some effects even when no specific states are chosen. 



4.3.4 Cosmic forgetfulness and signature change 

There are additional and more- surprising properties of the high-density regime that inval- 
idate a traditional bounce interpretation even in the harmonic model. First, staying in 
the isotropic context, there is cosmic forgetfulness [5SJ EH]: When crossing the bounce in 
^-evolution, some moments change in ways so sensitive to initial values that the pre-bounce 
state cannot be recovered precisely from what may be known post bounce. Using solutions 
of moment equations [89] or considerations of semiclassical wave functions [90J, one can 
derive an inequality 

(AVV 



< ^ilM (94) 
" (AV/(V))^ 
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bounding the ratio of volume fluctuations at early and late times. For a state with fluctua- 
tions symmetric around the high-density regime near the minimum of (V)(<f)), the left-hand 
side would be near zero. The estimate can therefore be used to shed light on the question 
of how much a quantum state may change while evolving through high density, and how 
much possible changes can be controlled. 

With matter fluctuations Ap^/ (p^) usually much larger than geometry fluctuations 
(AV/ (V))^oo at large volume, where quantum field theory on curved space-time should 
be a good approximation, the right-hand side of the inequality is much larger than one, 
and (AV^^-oo can differ significantly from (Ay)^. The inequality can be saturated 
by highly squeezed dynamical coherent states [89], showing that control on pre-bounce 
fluctuations cannot be improved unless states are restricted further, more strongly than by 
semiclassicality. 

Several classes of specific states, especially ones with weak correlations of the canonical 
variables — the volume and the Hubble parameter — show more-symmetric behavior of 
volume fluctuations. Most wave functions that can be constructed explicitly, for instance 
using sl(2, R)-coherent states based on the algebra (!85|) as used in [91], are only weakly 
correlated and do not show all possible asymmetries. Again, the effective viewpoint using 
moments instead of wave functions provides larger generality. And even though wave 
functions are not provided in explicit terms, one can show that wave functions for the 
moment solutions even at saturation of (194"1) do indeed exist: examples for such states are 
dynamical coherent states [57] saturating the uncertainty relation, whose existence can be 
shown by general methods well-known from quantum mechanics. 

The second feature preventing a bounce interpretation brings us back to quantum space- 
time structure. For a reliable cosmological model, we must embed a holonomy- modified 
isotropic version within a consistent deformation of the constraint algebra. Only then 
can we be sure that the model describes consistent evolution of quantum space-time. No 
complete extension of holonomy corrections to inhomogeneity is known, but as we will 
see in the next section, existing versions of holonomy modifications at high density imply 
drastic modifications with signature change, turning space-time into a quantum version 
of 4-dimensional Euclidean space. This happens right where the bounce would be, but 
without time and evolution in Euclidean space, a bounce interpretation is not valid even 
though the model remains non-singular. 

5 Quantum geometry and dynamics of space-time 

So far, in Section [3j we have seen the quantum geometry of space, with its characteristic 
features of discrete structures in loop quantum gravity. To fit these modifications into 
a covariant quantum space-time structure, we must find a consistent deformation of the 
hypersurface-deformation algebra of which the modified Hamiltonian constraint is a part. 
As has by now become clear from many examples, derived at different levels of effective and 
operator calculations, the classical constraint algebra is then indeed deformed: quantum- 
geometry corrections imply modified quantum space-time structures [llj. The possibility 
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of such deformations has also been suggested based on Wheeler-DeWitt quantization [92] . 
Instead of the classical commutator of two time deformations, we have 

{HimiHiNz]} = DlPq^NiVbNz - AW^i)] (95) 

with a phase-space function (3^1. 

The typical and rather universal form of these deformations is as follows. Holonomy 
corrections result in a curvature- dependent (3(K) = cos(2£K), where £, related to the 
quantum-gravity scale L, is a holonomy parameter depending on the curves used to inte- 
grate the connection, and K is a curvature component such as the Hubble parameter for 
perturbations around isotropic models [93J or the rate of change of orbit areas in spheri- 
cal symmetry [HU [95] . Such a deformation has also been found by operator calculations 
in 2 + 1-models |96j. For inverse-triad corrections, (3 = a 2 depends on the inverse-triad 
correction function a as in ( 1391) . which in turn depends on the quantum-gravity scale L 
[TT [ |94" 1 [95]. Also here, operator calculations have provided supporting evidence [55] . 



5.1 Example: spherical symmetry 

Spherically symmetric models with their reduced number of free fields provide an interest- 
ing testing ground for different quantum space-time structures, and at the same time allow 
physical applications for instance to black-hole physics. In Ashtekar-Barbero variables, 
used to compute deformations with corrections from loop quantum gravity, we express the 
canonical structure by four fields, two scalars A v and E x and two densitized scalars A x 
and E v . They appear as components of spherically symmetric SU(2)-connections 

A = A x (x)r 3 dx + A^(x)A A di!) + A (p (x)A A sin + r 3 cos#d<£ (96) 

and densitized triads 

E = E x (x)t 3 sin — + E*(x)A E sin tf — + E^(x)A E — . (97) 
ox ov dip 

For the general derivation see [97J [981 H] • 

The su(2)-valued fields A E = r!Cos(C(x)) + T 2 sm(((x)) and A E = r^A^ describe a 
U(l)-gauge freedom with gauge rotations by exp(A(x)r 3 ), remnant from the initial SU(2)- 
freedom. Similarly, the densitized triad has independent su(2)-matrices and A^. Also 
A x is affected by these gauge transformations, under which it changes to A x + dA/da; like a 
U(l)-connection, but the combination A x + d(/dx is invariant, and happens to agree with 
an extrinsic-curvature component K x (up to a factor of 7). 

With different matrices and A E , the components E v and A v , unlike E x and A x or 
K x , are not canonically conjugate. However, if we switch to extrinsic curvature also for ip- 
components, we obtain canonical pairs {K x (xi), E x (x 2 )} = 2'jG5(xi, x 2 ) and {^(^i), -E </3 (x 2 )} 
, yG5(xi,X2), as shown in [99]. In these canonical variables, we have the diffeomorphism 
constraint 

D grav [N x ] = j dxN x (2K'^ - K x E xl ) (98) 
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and the Hamiltonian constraint 

H[N] = ~^J dxN\ET 1/2 ((1 - + K%)E* + 2\E X \K^K X + 2\E X \T' V ) . (99) 

with = —(E x )'/2E tp a spin-connection component. (Primes denote derivatives by x.) 

The appearance of inverses of E x and quadratic expressions in K x and K v , the latter 
as canonical and gauge-invariant versions of the connection, suggests inverse-triad and 
holonomy corrections from loop quantum gravity. Hamiltonian constraint operators have 
been constructed in [HH], and more generally for Gowdy models in [100[ 1101] . in which 
inverse-triad operators and holonomy operators indeed appear. An effective Hamiltonian 
would be obtained from expectation values of these operators, but the calculations are 
complicated. Moreover, so far these Hamiltonians could not be ensured to be anomaly-free 
at the operator level. Instead, we can parameterize effective Hamiltonians by correction 
functions originating from inverse-triad and holonomy operators, and compute Poisson 
brackets of the modified constraints to see under which conditions they can be anomaly- 
free jgu. 

We write a general modified Hamiltonian constraint as 

H%JN] = -^JdxN( y a\ET 1/2 E ifi f 1 (K ip ,K x ) + 2a\E x \ 1 / 2 f 2 (K^K x ) 

+a r |£T 1/2 (l - K) EV + 2 ®r^\E x \ 1/2 ^ , (100) 

with inverse-triad correction functions a, a, ar and ar initially left independent of one 
another, and holonomy correction functions f\ and f 2 . All these functions may in principle 
depend on all canonical variables, although the triad dependence of inverse-triad correction 
functions and the curvature dependence of holonomy correction functions should be pri- 
mary. Moreover, an anomaly-free commutator with the diffeomorphism constraint shows 
that inverse-triad correction functions can only depend on E x , not on the density- weighted 
E*. 

Computing Poisson brackets of two modified Hamiltonian constraints with different 
lapse functions, it turns out that anomaly freedom can be realized if f\ = F 2 and f 2 = K X F2 
provided that F 2 = Fi(dFi/dK v )a/ar [94J. Choosing a function Fi periodic in K v , holon- 
omy modifications for this component are realized. The second correction function F 2 
is then fixed, showing how anomaly-freedom can put restrictions on possible modifica- 
tions and quantum corrections. In fact, the corrections seem even stronger for the in- 
dependence, left unmodified in the function f 2 shown here. An extension to a holonomy- 
corrected ^-dependence appears more difficult than one of the ^-dependence. Moreover, 
while K v would give rise only to pointwise exponentials exp^iK^) as holonomies, the 
curve integration along angular directions, in which K v points, being trivial in spherically 
symmetric models, K x would be replaced by a holonomy exp(«7 fjdxKx) integrated along 
some interval /. Derivative corrections should therefore result as well, or the constraint 
would become non-local if integrations are left unexpanded. 
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If we take Fi(K v ) = (^i)" 1 sm^jiKp), suitable for holonomy corrections as in the cos- 
mological example (166|) . we have F 2 (K ip , E x ) = (2 r y£)~ 1 sm^^iK^a/ar- The algebraic 
deformation is then given by (3(E X , K v ) = aardF 2 /dK lfi [94J. For the example provided, 
this means f3(E x ,K ip ) = a(E x )a(E x )(ctr(E x ) / cer(E x )) cos(2'j£K ip ), a function that is neg- 
ative for 'jiKp ~ 7r/2, at curvatures where the correction function fx is near its maximum 
and a strong modification of the classical linear function. This property is realized generi- 
cally: Combining the previous equations, we can write 



which is negative around maxima of fx, irrespective of its functional form. Consequences 
of negative f3 will be discussed in more detail soon. Note also that holonomy corrections 
and inverse-triad corrections are rather independent of each other in their effect on the 
deformation, affecting /3 multiplicatively. 

The inverse-triad correction functions must satisfy 



for a closed constraint algebra [94]. If F\ is independent of E x , or at least depends on this 
triad variable in a way different from inverse-triad corrections, one can show that both 
terms in (11021) must vanish individually, and we have «r = ot and «r = a. See also [102J. 

For consistent deformations in the presence of cosmological perturbations, anomaly- 
freedom is implemented in the same spirit, but with an extra ingredient. Without any 
symmetry assumptions, requiring irregular lattices and non-Abelian SU(2)-features, it is 
complicated to derive inverse-triad operators or to parameterize holonomy corrections. 
One therefore starts using all information about such correction functions that can be 
obtained in tractable models, such as a homogeneous background, and inserts those back- 
ground functions just like a, a, ar, ar, fx and f% in spherically symmetric models. These 
functions refer only to the background variables, but depending on the order of cosmolog- 
ical perturbations, also the dependence on inhomogeneity is required. The corresponding 
terms, in many cases, cannot be computed directly from operators; instead, one inserts 
"counterterms" in the Hamiltonian constraint expanded by inhomogeneity, taking into ac- 
count all possible terms to the given order that could be generated by correction functions 
depending on homogeneous fields [11]. Terms that cannot be computed from operators 
are left unspecified as free functions. In many cases, counterterms contribute derivative 
corrections, adding for instance terms containing d a E\ for inverse-triad corrections. A de- 
pendence of correction functions on integrated variables such as fluxes is therefore realized 
even if a initially depends only on local triad values. The condition of anomaly-freedom is 
often so restrictive that the counterterms can be derived uniquely from known inverse-triad 
or holonomy correction functions of the background. Consistent constraints can therefore 
be computed even if not all quantum corrections are known in detail. 

Another small difference between spherical symmetry and cosmological perturbations 
is the treatment of the SU(2)-gauge. In spherical symmetry, ( 198|) and fllOOp are manifestly 




(101) 




(102) 
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invariant under these transformations, even after deformation by quantum corrections. In 
perturbative treatments, on the other hand, one usually fixes a background triad, including 
its SU(2)-gauge. But then, one can see that physical results and the deformation of the 
space-time algebra are independent of the specific SU(2)-fixing chosen. The constructions 
are therefore consistent even if some gauge has been fixed: All gauge fixings produce the 
same results. (The same arguments can be applied to the time gauge used to descend 
from space-time tetrads to spatial triads.) Such consistent derivations in the presence 
of gauge fixing are possible in simple cases such as the Gauss constraint, which moreover 
survives unmodified after quantization. Making gauge-fixings consistent is much less trivial 
if attempted for complicated constraints such as H[N}. First, it is difficult to find any 
good gauge fixing in general terms; having to study even all possible gauge fixings and to 
make sure that physical results do not depend on the choice is then nearly impossible. In 
such cases, the only manageable approach is to forgo gauge fixing before quantization or 
deformation, even if it makes derivations more complicated than in one given gauge. 

5.2 The meaning of deformed hyper surface deformations 

The hypersurface-deformation algebra encodes the space-time structure of generally co- 
variant theories just as the Poincare algebra encodes special relativity's structure. One 
can recover the Poincare relations by using functions N and N a linear in some coordinates 
amounting to Minkowski space-time or a local Minkowski patch. For instance, two linear 
lapse functions of the form N = At + v ■ Ax, inserted in (1121) . provide the commutator 
of Lorentz boosts by velocity v and time translations At. Inserting two such functions, 
Ni = v- Ax and N2 = At — v ■ Ax, on the right-hand side of (fl2l) shows a commutator that 
amounts to the displacement Ax = vAt. In terms of linear hypersurface deformations, the 
relation follows from elementary geometry; see Fig. [TJ 

If the algebra is deformed, as in ffl3|) . the same choice of linear deformations along the 
normals gives rise to a rescaled relation Ax = (3vAt. Quantum space-time, with its discrete 
structure that is responsible for the algebraic deformation via holonomy and inverse-triad 
corrections, changes the relation between velocity and displacement; discrete space-time 
speeds up or slows down motion. Such a phenomenon is well-known from condensed-matter 
physics and should not come as a surprise, although the form in which it is realized here 
is rather different owing to the more-basic notions of space and time involved. 

When (3 becomes negative, as happens with holonomy corrections at high density, 
the relation Ax = (3vAt is rendered counter-intuitive. However, the change of sign can 
be interpreted easily if one redraws Fig. [T] in Euclidean space, especially regarding the 
directions of normals to spatial slices. As shown in Fig. [2] compared with Fig. [TJ the 
displacement then indeed points in the opposite direction (and does not change magnitude). 
The relation Ax = — v At does not describe motion but rather, despite the notation of 
variables, a rotation. When /3 turns negative, we have Euclidean signature rather than 
Lorentzian. Even if /3 is not exactly —1 but negative, the structure is best described as 
Euclidean even though we do not have classical Euclidean space (just as we do not have 
classical Minkowski space if (3 is positive but not exactly +1). 
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Figure 1: Two linear deformations of spatial slices, first by N\ = vAx along the normals 
and then by N 2 = At — vAx (top) and in the opposite ordering (bottom), commute up to a 
spatial displacement by Ax = vAt. Normals are drawn according to Minkowski geometry, 
corresponding to Lorentzian signature. 




Figure 2: With normals drawn according to Euclidean geometry, the displacement seen 
in Fig. CD points in the opposite direction: The sign in the commutator of two normal 
deformations is reversed. 
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Signature change is not only a drastic reminder that we cannot take much of our 
usual concepts for granted when space-time is quantized. It also shows the limitations 
of approaches in which gauge-fixing or deparameterization is used before quantization or 
modifications of the constraint. If one fixes the gauge before quantization, one can only 
assume the classical space-time structure, which does not allow signature change. After 
gauge fixing and quantization, the full space-time structure can no longer be accessed, leav- 
ing one with the conclusion that space-time is unmodified. However, since the constraints 
are modified in quantization after gauge fixing and determine the gauge and space-time 
structure, the procedure becomes intrinsically inconsistent. Deparameterization cannot 
capture all quantum space-time effects either. When one distinguishes a phase-space de- 
gree of freedom to measure change and evolution, there is no guarantee that this degree 
of freedom actually behaves like time in a space-time sense. Also Euclidean theories can 
formally be deparameterized (internal time simply parameterizes gauge orbits of the Hamil- 
tonian constraint), showing that deparameterized "evolution" does not necessarily imply 
evolution in a temporal sense. Again, only a complete analysis of the off-shell constraint 
algebra, without eliminating some of its more complicated ingredients by gauge-fixing or 
deparameterization, can show what space-time structure is realized. 

In our diagrams so far, we have assumed that normals are drawn using either Minkowski 
or Euclidean geometry. We did not use quantum corrections of angles even though dis- 
tances and displacements did receive corrections as a consequence of the algebra ( fl3l) . One 
may expect angles to change too, in particular angles of normals to spatial slices drawn to 
visualize the commutator of two time deformations. Such corrections could indeed happen, 
in general deformations of space-time structures, but since the angle between space and 
time directions would be involved, they would amount to deformations of the commutator 
( II ip of a time and a space deformation. A temporal and a spatial deformation commute up 
to a temporal deformation, as illustrated in Fig. El If there are quantum corrections to this 
relation, one might interpret them as modified spatial displacements, rescaled compared 
to the classical relations, or a modification in space-time angles used to define temporal 
deformations along the normals. In the former case, also the purely spatial commuta- 
tor ( ITU]) should be modified, which is not the case. A modified (JTTJ) therefore indicates 
quantum corrections to space-time angles, not just to distances as indicated by a modified 
( Tl2|) . For instance, modifications which add a term of D[(3'NN a ] to the classical result of 
{H[N], D[N a }}, as possible for vector modes [103J but not with scalar modes [93], would 
have an additional spatial shift of the open circles in fl3]). (Note that Fig. [3] also illustrates 
the fact that (1111) is not subject to a change in sign in Euclidean signature. If the normals 
are drawn according to Euclidean geometry, just the rescaling of Ax under boosts changes, 
but not the direction of the temporal displacement.) 

Based on algebraic calculations in effective loop quantum gravity, modifications of (fTTj) 
that seemed possible for vector modes subject to holonomy corrections |103] turned out 
not to be consistent with scalar modes [93]. It therefore appears that space-time angles are 
not affected by the corrections of loop quantum gravity, but as discussed in Section 13.21 
stronger corrections to the diffeomorphism constraint than used so far may be realized. 
For now, all consistent deformations point to quantum corrections only in the commutator 
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Figure 3: The commutator of linear time and space deformations: A Lorentz boost pro- 
duces the solid inclined line, the same boost performed after a spatial displacement the 
dashed line. A spatial point (full circle) is mapped to two different positions (open circles) 
depending on whether the spatial displacement is performed before or after the boost. The 
final positions differ by a time translation (dotted), according to {H[vx}, D[Ax]} = H[vAx] 
using ( fTTj) . 



of two time deformations. 

5.3 Consistency and space-time structures 

With a deformed hypersurface-deformation algebra the theory is fully consistent. The 
full amount of gauge transformations exists to remove spurious degrees of freedom, and 
constraints are guaranteed to be preserved by evolution. Cosmological observables, for 
instance, can be computed by standard Hamiltonian means, as developed for deformed 
algebras in (TTJ 1104] . (The classical Hamiltonian formalism for cosmological perturbations, 
of |105j . assumed certain features of observables that are modified by deformations.) 

However, other familiar notions used often in general relativity no longer apply. Even 
the line element, one of the most basic mathematical objects of differential geometry, 
must be treated with care — or altogether avoided. Constraints obeying a modified al- 
gebra generate gauge transformations for metric components q a f, that cannot agree with 
Lie derivatives or space-time coordinate transformations — otherwise, the form of gauge 
transformations would imply the classical hypersurface-deformation algebra. But if qw,, or 
the space-time metric g^ v completed in the usual way using lapse and shift in 

ds 2 = g^dx^dx" = -N 2 dt 2 + q ab (dx a + N a dt)(dx b + N b dt) , (103) 

does not transform by classical coordinate changes, the contraction g^dx^dx" with stan- 
dard coordinate differentials dx M is not invariant and cannot be used as a line element. 
One would have to modify transformations of coordinate differentials as well to make the 
contraction invariant, in which way one could possibly make contact with non-commutative 
[106] or fractional geometry |107j . (Modified hypersurface-deformation algebras have also 
been found with higher-derivative dispersion relations for matter [108J, but the form does 
not appear related to what is suggested by loop quantum gravity.) 
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Lacking a line element and related notions such as geodesies or trapped surfaces, black- 
hole properties and even their definitions must be rethought. One can only rely on canonical 
formulations of horizon conditions that capture the well-known notions of standard space- 
time. Classically, there are several different definitions of horizons which, at least in simple 
cases, all provide the same results. However, when they are reformulated canonically and 
adapted to quantum space-time, results may differ. As analyzed in (109J, the consistency 
of quantum space-time structures also in this case helps to arrive at unambiguous answers. 

For calculations in quantum space-time, one must rely on Hamiltonian methods, using 
directly the modified constraints. First applied systematically in |104] , one obtains unam- 
biguous expressions for gauge-invariant variables in cosmology as well as their evolution 
equations. Without using an action, it may not always be obvious how to combine all 
equations obtained to just one Mukhanov-type equation for the analog of the curvature 
perturbation, but examples have been found for inverse-triad corrections [HOj and holon- 
omy corrections [111] . Quantum corrected wave equations are then obtained, which directly 
show the modified speeds of modes as expected from a deformed hypersurface-deformation 
algebra: electromagnetic and gravitational waves obey the equation 



with the correction function /3 from ( |T3|) and a function / that shows how the evolution of 
modes depends on an expanding background. Density perturbations obey a similar equa- 
tion, but with a differently modified speed for inverse-triad corrections [110J. Interesting 
phenomenological effects are therefore suggested. (For holonomy corrections in currently 
existing versions, the modified speeds of density perturbations and gravitational waves are 
identical [TT2] .) 

Holonomy modifications at high density are especially drastic: they imply signature 
change |102l 1113] . We are no longer dealing with quantum space-time but with a quantum 
version of 4-dimensional Euclidean space, shown by the sign change of 0(K) for large 
IK. The deformed hypersurface-deformation algebra (TLB"]) then belongs to Euclidean-type 
space, and (11041) shows the elliptic nature of linear mode equations. For positive as well 
as negative {3, we do not have standard space or space-time unless the value is exactly ±1; 
quantum space-time effects always occur. But the distinction between positive and negative 
values of (3, as opposed to two different positive values, is much more important because 
it changes the type of initial-value or boundary problems in mode equations. Signature 
change defined by the sign of j3 is therefore physically relevant, even in the absence of a 
standard classical space-time structure. In these high-density regimes, however, it is no 
longer consistent to treat holonomy corrections in isolation because they mix with quantum 
back-reaction, together forming higher-curvature corrections. 

So far, no consistent deformation of the constraint algebra has been found for quantum 
corrections caused by quantum back-reaction, adding moment-dependent terms to the 
classical constraints. The problem is well-defined because we know the Poisson algebra of 
moments, and even though canonical effective field theory techniques remain incomplete, 
one could use those of quantum mechanical models in loop quantum gravity restricted to 




(104) 
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fixed graphs. In any finite region, there would then be a finite, though large, number of 
degrees of freedom given by link holonomies and plaquette fluxes. (Most calculations in the 
full theory are done with fixed graphs, anyway.) Since Poisson brackets of moments always 
produce other moments, quantum corrections must be arranged so that surplus terms 
delicately cancel in the constraint algebra. No such version has been found yet, which is 
not altogether surprising because a successful implementation would imply a consistent 
version of quantum space-time, including moments of a state preserving covariance. 

5.4 Anomaly problem 

The problem of finding consistent deformations of the hypersurface-deformation algebra 
is related to the long-standing anomaly issue of canonical quantum gravity. If one can 
find an anomaly-free representation of quantum constraint operators, turning the classical 
constraint algebra into some operator version, effective constraints will automatically be 
first class and consistent: If two constraint operators C\ and C 2 commute up to another 
constraint operator, the quantum constraints (Ci) and (C 2 ) Poisson commute up to another 
quantum constraint, thanks to {(Ci), (C 2 )} = ([Ci, C 2 ])/ih. This statement also extends 
to higher-order quantum constraints (f(q,p)C) [37] . 

Structure functions in the constraint algebra cause several problems in trying to find 
anomaly-free operator versions. But if such a version has been found, there are no addi- 
tional problems in the transition to effective constraints. The product [C\, C 2 ] = fC with a 
quantized structure function / will, after taking an expectation value, be one of the higher- 
order quantum constraints. However, there are obstructions to simple attempts at solving 
the anomaly problem for constraint operators, related for instance to the interrelation of 
ordering and self-adjointness issues especially in the presence of structure functions |114j . If 
Hamiltonian constraint operators are self-adjoint, the commutator [H[M], H[N]]/ih quan- 
tizing {H[M}, H[N]} in f|T2|) is self-adjoint too. But classically the bracket equals a product 
of the local diffeomorphism constraint D a with the structure function q ab , two expressions 
that have a non-vanishing Poisson bracket because 

q ab ig 

not diffeomorphism invariant. 

Any quantization of the product q ab D a can then be self-adjoint and have a chance of 
agreeing with [H[M], H[N]]/ih only if a symmetric ordering is used. But if one simply 
reorders u q ab D a " symmetrically, a metric factor would come to lie to the right of the dif- 
feomorphism constraint, and the product would no longer annihilate physical states. An 
anomalous version of the constraints would be obtained. 

In an effective description, the situation is more manageable. First, one can use cal- 
culations of Poisson brackets instead of commutators of operators, even in the presence of 
quantum corrections and ordering choices. Moreover, quantum corrections may be suitably 
parameterized, to take into account different ordering choices, quantization ambiguities in 
the representation of inverse-triad and holonomy operators, and general classes of states in 
terms of their moments. One can compute Poisson brackets with ambiguity functions, such 
as (159]) . or moment terms unspecified, and see what conditions anomaly-freedom imposes 
on them. An operator calculation with free functions or states, by comparison, would be 
much more involved. The self-adjointness question can be left open at first by allowing 
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for complex-valued constraints — kinematical moments appearing in effective constraints 
are complex-valued anyway, even in the absence of structure-function issues, as seen in 
Section 14.1.51 With this strategy the first consistent deformations of the form (fT3l) have 
been found by now confirmed also by operator calculations [55]. The effective view is 
therefore reliable and powerful also in the context of the anomaly problem. 

5.5 Quantum back-reaction and higher time derivatives 

Moments and their quantum back-reaction, if they appear in a consistent deformation of 
the constraint algebra, are the canonical analog of higher-curvature terms with their higher 
time derivatives. However, their form is not directly one of higher time derivatives, and 
additional steps and expansions, primarily an adiabatic one, are required to put canonical 
effective equations in the form of higher-derivative ones; see Section 14.1.21 The adiabatic 
approximation is not always applicable. It serves well for anharmonic oscillators expanded 
around the harmonic vacuum with its constant moments. Moments of the anharmonic 
vacuum change, but only slowly, and can be treated adiabatically. Solving for moments 
order by order in the adiabatic expansion then shows how they are related to higher time 
derivatives of expectation values. The same statements hold true in quantum field theory, 
where one expands around the free vacuum in order to describe excitations around the 
interacting vacuum using the low-energy effective action. 

These features have given rise to the expectation that quantum gravity should be sub- 
ject only to higher- curvature corrections. However, two hidden assumptions are required 
for this conclusion. First, one assumes that quantum gravity implies corrections only in the 
dynamics of gravitons, say, not in the underlying space-time structure. Since gravity theo- 
ries are fundamentally about space-time structure, this assumption, valid in a perturbative 
context of excitations on a fixed background, need not be true in general. 

The second assumption is that quantum gravity can be realized perturbatively around 
some free theory of the usual form. Also this statement may be true for perturbations 
around a fixed background and perhaps some other situations, but it is not always valid. 
Quantum gravity is not expected to have a non-perturbative ground state, and other 
distinguished states may be very different from Gaussians as they appear in the harmonic- 
oscillator ground state or free vacuum states. Such states may not allow adiabaticity of the 
moments, and therefore cannot have a dynamics fully expressed by higher time derivatives. 
Indeed, examples in quantum cosmology are known in which the adiabaticity assumption 
is difficult to realize consistently |115j . Even for anharmonic oscillators, the adiabatic 
approximation is not valid when one perturbs around correlated coherent states, such as 
fully squeezed Gaussians, of the harmonic oscillator instead of the ground state [301 170] . 

Instead of a higher-curvature or other higher-derivative effective theory, canonical quan- 
tum gravity has higher- dimensional effective systems in which the moments play the role 
of independent degrees of freedom in addition to expectation values. They are subject 
to their own equations of motion, and by quantum back-reaction couple to them, as e.g. 
in ( H6|) - (l50|) . Their effect can be formulated by higher-derivative terms only in certain 
regimes, but not in general. 
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5.6 Effective actions 



So far, we have stayed in the canonical framework and dealt with effective equations and 
constraints. Complementary information can be obtained if a corresponding effective ac- 
tion is found. In principle, there is a one-to-one correspondence between canonical and 
Lagrangian formulations, but in perturbative settings, especially those that imply higher 
derivatives, the transformation is far from obvious. A derivative expansion of a Lagrangian 
may then correspond to a complicated resummation of a derivative-expanded Hamiltonian. 
It is therefore of interest to look for and study effective actions even if effective canonical 
equations are already known. 

There are several examples in which effective actions have proven their usefulness in 
quantum cosmology. Euclidean path-integral techniques have been developed and ap- 
plied to questions such as tunneling probabilities and some semiclassical issues |116[ I117[ 
11181 1119] . As expected, in semiclassical regimes, corrections to Einstein's equation are of 
higher-curvature form. Also causal dynamical triangulations have led to effective actions of 
cosmological systems by comparing detailed numerical studies of volume fluctuations with 
the dynamics of minisuperspace models [120] . By matching volume correlation functions 
with results expected from a higher-curvature effective action, quantum corrections can 
be derived. A general comparison involving also possible modifications to quantum space- 
time structure on top of higher-curvature corrections has not yet been completed. Effective 
actions of such forms often show more directly than other types of effective equations when 
quantum effects become significant. 

Effective actions based on path integrals are subject to modified quantum space-time 
structures just as canonical effective equations and constraints. Path integration for grav- 
itational theories requires an integration over all metrics, with a suitable measure that 
preserves covariance and does not introduce anomalies. Such a measure has not been 
found in complete generality, indicating that one encounters in this approach the same 
difficulties that appear when one tries to represent the canonical constraint algebra con- 
sistently. The path-integral measure is indeed the place to look at when one is interested 
in possible space-time modifications, but its incomplete nature does not give many clues. 
Canonical theories, with consistent deformations of the constraint algebra found in recent 
years, have been able to make more progress in this direction. (The spin-foam approach 
|121[ I122[ 1123] attempts to take the results of loop quantum gravity regarding background- 
independent representations to a level comparable to path-integrals [124] . However, while 
some aspects of integrations over the space of metrics can be clarified, the issue of the cor- 
rect measure, or spin- foam face amplitudes, remains open also here |125j . Notwithstanding 
these problems, traditional effective-action techniques have been applied to spin foams in 
[126].) 

Another question in path-integral based approaches is how the state dependence enters 
effective actions. The low-energy effective action can be computed quite conveniently with 
path integrals, but it hides the fact that one is expanding around the ground state. If one 
tries to go beyond this limitation, which in quantum cosmology with its lack of ground 
states is a severe one, the required calculations become much more involved. Suitable states 
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would have to be implemented by additional wave-function factors in the path integral, and 
integrations would no longer be Gaussian for general states. Moreover, since path integrals 
in their most common form provide transition amplitudes between two states separated 
by some finite time, one would have to put in the initial and final state. Compared to 
effective canonical equations, which only require initial moments to be specified, more a- 
priori information about the physical system is therefore required. In quantum cosmology, 
this information is not easily found by independent means, but is important for an analysis 
especially of the Planckian regime. 

5.7 Regained dynamics: from canonical to Lagrangian 

In this situation, a combination of canonical and Lagrangian effective methods is of inter- 
est. As already mentioned, it may be difficult to perform a transformation for constrained 
systems, especially those with higher-derivative terms. If the relation between the momenta 
and time derivatives of configuration variables involves higher derivatives, as easily happens 
with quantum-gravity corrections such as holonomy modifications, the Lagrangian would 
appear as a complicated resummation of the higher-derivative expansion of the Hamilto- 
nian. Moreover, even if one can Legendre transform the Hamiltonian constraint to arrive 
at a Lagrangian, one would, for a comparison with path integrals, still have to look for 
a measure under which the modified Lagrangian is covariant. Also integrating the La- 
grangian to an action requires new constructions even though just a classical measure on 
space-time is needed. However, with space-time structures modified and effective metrics 
non-existent — see Section 1531 — defining space-time integrations is non-trivial. Measure 
issues, therefore, cannot be resolved easily, but at least an effective Lagrangian for further 
classical-type analysis may still be found. 

Instead of attempting a Legendre transformation, the canonical methods of [127} 1128] 
can be used to derive a Lagrangian (or the constraints themselves) directly from a modified 
constraint algebra, to any given order in derivatives. If only second-order derivatives are 
assumed and the constraint algebra is classical, the Einstein-Hilbert action as a two- 
parameter family with Newton's and the cosmological constant is obtained as the unique 
solution. If higher derivative orders are allowed, one expects higher-curvature corrections 
as well. And if even the constraint algebra is quantum corrected, as suggested by loop 
quantum gravity, stronger corrections, also to second derivative order, are obtained. With 
this procedure of "regaining" a Lagrangian from the constraint algebra one can sidestep 
the complicated resummations that a Legendre transformation from modified constraints 
to the higher-derivative Lagrangian would imply. Not all coefficients in the Lagrangian 
may follow uniquely, especially if higher derivatives are included leaving several options of 
higher-curvature invariants of the same order. But the general form of modifications and 
implications for quantum space-time structure can still be found. 
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5.7.1 Functional equations 



To spell out the general procedure, let us assume that we have a Hamiltonian or Hamilto- 
nian constraint H(q,p) depending on canonical fields q(x) and p(x) in space. We introduce 
SH/8p(x) =: v(x) as a new independent variable in place of p, and then expand equations 
by this newly defined v. If H is a Hamiltonian generating evolution in some fixed time 
parameter, we have v(x) = q(x) by Hamilton's equations. If if is a Hamiltonian con- 
straint in the absence of an absolute time, v = (5N)~ 1 {q, H[5N]} is the derivative of q in 
a direction normal to spatial slices. For gravitational variables with q the metric or triad, 
v would therefore be related to extrinsic curvature. This change of variables amounts to 
what is needed for a Legendre transformation from (q,p) with Hamiltonian H to (q, v) with 
Lagrangian L = pv — H, whose form will result as a solution of the f-expanded equations. 

Note that we cannot always assume the Hamiltonian to be local and free of derivatives 
of p, which would imply that partial derivatives could be used to compute v. Holonomy 
corrections in loop quantum gravity, for instance, introduce higher spatial derivatives of 
the momentum conjugate to the densitized triad. In such situations, there is no local 
relation between the Hamiltonian and the Lagrangian, and explicitly performing a Legendre 
transformation is complicated. 

Instead of computing the transformation, we intend to calculate the Lagrangian di- 
rectly from the constraint algebra, assuming from now on the relations of the (deformed) 
hyp ersurface- deformation algebra. Using the definition of v and 
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Sq(x') 



_ SL 



(105) 

v(x) 



for the unsmeared Hamiltonian constraint and the Lagrangian density, we write the Poisson 
bracket ( !T3|) of two Hamiltonian constraints as 

{H[N],H[M]} = - f d 3 x f d 3 y^^-v(x)N(y)M(x) - (N M) (106) 

J J oqyx) 

d 3 x(3D a (x)(NV a M - MV a N) (107) 



with the local diffeomorphism constraint D a . Taking functional derivatives by iV and M, 
we arrive at the functional equation 

^f^-v(x') + f3{x)D a {x)V a 5{x, x 1 ) -{x^x') = (108) 
oq[x ) 

for L(x), which can be solved once an expression for the diffeomorphism constraint D a is 
inserted, depending on whether q refers to gravity or some matter field. In all cases, D a is 
linear in the momenta. A linear equation for L is thus obtained |128j . If fflOl) and (fTTl) are 
unmodified, standard expressions for D a can be used. (See Section [32] for a discussion of 
possible further modifications.) 
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5.7.2 Matter 



To illustrate the regaining procedure, we look at a scalar matter field 4> without derivative 
couplings, whose Hamiltonian obeys the (deformed) hypersurface-deformation algebra on 
its own, without adding the gravitational piece. The Lagrangian density must be of the 
form C = y/det gL(<f), v, ip) where v = (5N)~ 1 {(p, H[5N]}, as before, is the normal scalar 
velocity and tp = g ab V a 0Vb0 is the only remaining scalar that can be formed from <f) and 
its derivatives, to a total derivative order of at most two. Higher derivatives may easily 
result in interacting matter or quantum-gravity theories, with higher time derivatives from 
quantum back-reaction and higher spatial ones, additionally, from possible discretizations. 
(See e.g. |129[ I130[ I131[ I132] for information about discretized space-time theories.) For 
now, however, we look for modifications implied by corrections that leave the classical 
derivative order unchanged. 

With the canonical variables of a scalar field and its diffeomorphism constraint D a = 
p^V a 0, Eq. ( I108p assumes the form 

5 -^TtMx') + /3^M(V>(x))V a 5(x, x>) - (x x') = . (109) 

As in [128] . we write 

SL(xl dL{xlH{x)_ 2 M£) 

5<j>(?) dcj>{x) 8cf>{x') + dij(x) [V WJW.^- 

It follows that 

A* := ( VV) (P^ + 2v^f] 



dv dif) J 

satisfies the equation A a (x)'V a 8(x, x') — (x 4-> x') = 0, shown in |128j to imply A a = 0. 
Thus, 

dL +2y dL _ Q 
dv dip 

and L must be of the form L(<p, ip — v 2 / (3). With non-trivial deformation, /3 ^ 1, the scalar 
field therefore obeys a modified dispersion relation. The kinetic term of the Lagrangian does 
not depend on ip — v 2 = (? Miy (V At 0)(V y 0) in space-time terms, but has its time derivatives in 
ip — v 2 / j3 rescaled by the correction function (3. The resulting modified dispersion relation 
is in agreement with the wave equation (11041) . 

At the canonical level, for comparison, we begin with a matter Hamiltonian density of 
the form 

H = v - P ±— + 1-ay/d^tP + y/feX~q W((f>) (110) 
zydet q 2 

with general inverse-triad correction functions v and a, and some potential W(<f>). The 
corresponding Lagrangian density, with v = up^/y/det q, takes the form 

L = Vd^(^-^-^)j =- > /drti|^-^) - V^qW{<j>). (Ill) 
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This function has the same kinetic dependence as derived above, provided that (3 = z/cr, 
exactly the requirement for an anomaly-free constraint algebra in the presence of inverse- 
triad corrections, where /3 = a 2 [XT] . Purely canonically, this consistency condition can be 
seen to ensure causality in the sense that gravitational waves on quantum space-time travel 
at the speed of light, modified by a factor of \f\j3\ compared to the classical speed |133j . 
No super-luminal propagation happens when anomaly-freedom is taken into account, even 
if propagation speeds may be larger than the classical speed of light for a > 1. 

5.7.3 Effective action for inverse-triad corrections 

For gravity, one example has been worked out in quite some detail with these methods: 
inverse-triad corrections of loop quantum gravity. These corrections have a characteristic 
component independent of higher derivatives, and therefore can be analyzed already at the 
level of second-order equations. The resulting second-order effective Lagrangian, regained 
from a modified constraint algebra with a correction function independent of curvature 
components, is 

with "velocities" v a b, defined again as normal derivatives [102] . For classical gravity, v a b = 
2K a b would be proportional to extrinsic curvature. Compared with the classical action 
obtained for /3 = 1, the notion of covariance has changed: Space and time derivatives 
are corrected by different coefficients of ( 3 )-R and |/3| _1//2 sgn(/3) of \{y a bV ab — v®v\), 

respectively. This result is consistent with the fact that the underlying constraint algebra 
is modified, taking the form f |13|) . For this reason, as already mentioned in the context of 
line elements, the manifold structure required to integrate Lp to an effective action is not 
clear. 

Nevertheless, the effective Lagrangian shows several characteristic effects. First, inverse- 
triad corrections, having implications even without higher-derivative or higher-order terms 
in v a b, can easily be separated from holonomy effects and higher-curvature corrections. 
They are especially significant for small fluxes, where /3 becomes small. With these con- 
structions, it is possible to use inverse-triad corrections even when they imply strong 
modifications, regimes in which one would otherwise have to include full holonomy and 
higher-curvature corrections as well. These latter types of corrections are indeed present, 
but affect only higher orders in the f-expansion. Inverse-triad effects up to quadratic order 
in v are reliable even if the other corrections are not known precisely. 

Time derivatives are then dominant compared to the spatial Ricci scalar when /3 ap- 
proaches zero at small fluxes, indicating that near-singular geometries are controlled by 
homogeneous dynamics, strengthening the classical BKL scenario |134j by a no- singularity 
scenario in loop quantum gravity |135j. While holonomy effects cannot easily be seen at 
this level since they would manifest themselves only at higher orders in v a b and mix with 
higher-derivative terms, they have one drastic implication at high density if they are dom- 
inant. Then, the holonomy correction function j3 becomes negative, a feature taken into 
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account by the sign factors in (I112p . When the sign changes, the signature does too: at 
the Lagrangian level, one can interpret the effect by turning t into it in (11121) with v a t, in- 
terpreted as first-order time derivatives. (The constraint algebra again provides a rigorous 
interpretation of this signature change; see Section l5\2l ) At Planckian densities, holonomy 
effects are so strong that they turn space-time into a quantum version of 4-dimensional 
Euclidean space, lacking time and evolution. Temporal interpretations of high-density 
holonomy implications, such as bounces, are incorrect. Only large higher- derivative terms 
could prevent /3 from turning negative, but then other holonomy effects such as bounces 
would go away too. 

6 Implications for phenomenology and potential tests 

It is difficult to test quantum cosmology by any observational means, and given the sub- 
stantial lack of control over deep quantum regimes, devising high-density scenarios of the 
universe remains a highly speculative exercise. Higher-curvature corrections, as they always 
appear in quantum gravity and cosmology except in the most simple harmonic models, are 
not relevant at currently accessible scales. They are certainly important in the Planckian 
regime, but then the present theories are so uncontrolled that it is impossible to derive 
clear effects, and even if some could be suggested, they would most likely be washed away 
by the immense amount of subsequent cosmic expansion. 

If one goes beyond higher-curvature corrections as computed for Wheeler-DeWitt quan- 
tum cosmology in |136j . using for instance quantum-geometry effects from loop quantum 
gravity, the situation has a chance of being more optimistic. Several investigations have 
been performed pH Q3S1 EM EM HHl Q321 CGEl Q3SQ35], but not all are based on con- 
sistent implementations of inhomogeneity. The general picture is therefore still incomplete. 
One type of corrections, holonomy modifications, provides contributions of higher powers 
of the connection or extrinsic curvature, very similar to some parts of higher-curvature 
corrections. Holonomy corrections therefore cannot be separated from higher-curvature 
effects, and cannot provide more-sizeable consequences regarding observations. 

Moreover, the quantum space-time structure corresponding to holonomy corrections 
remains incompletely understood. Consistent deformations of the classical constraint al- 
gebra with some holonomy-like effects are known in spherically symmetric models [M| [95], 
2 + 1-dimensional gravity [96] and for cosmological perturbations [93J . However, in spher- 
ically symmetric and cosmological models, only "pointwise" holonomy modifications have 
been implemented, replacing connection components c by exp(i£c) but not integrating 
over curves. Curve integrations, on the other hand, provide additional terms which are 
non-local or, when a derivative expansion (|34|) is used, introduce higher spatial deriva- 
tives. Keeping only higher powers of c as in an expansion of the exponential, but ignoring 
spatial derivatives is not a consistent approximation: In strong curvature regimes, where 
holonomy effects should be significant, higher powers and higher derivatives of connection 
components both contribute to the same order of curvature. 

Inverse-triad corrections, fortunately, are much better-behaved. First, the derivation 
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of consistent deformations based on them is more complete, achieved in the same kind 
of models - - cosmological perturbations [11], spherical symmetry [94"l [95] . and 2 + 1- 
dimensional models [55] in which operator calculations can be performed - - but with 
all crucial terms included. Also inverse-triad corrections should come along with higher- 
derivative terms because they depend on fluxes, or integrated densitized triads. However, 
the non- derivative contribution is significant, as seen in Section l5.7.31 and, unlike derivative 
terms, depends on parameters unrelated to curvature components. Moreover, additional 
derivative terms are included by the counterterms of [llj, while no connection-derivative 
counterterms have been used for holonomy corrections in [93] . Inverse-triad corrections are 
therefore more reliable than holonomy corrections at the present stage of developments in 
loop quantum cosmology. 

Inverse-triad corrections are also more interesting from an observational perspective. 
Because they do not directly refer to the curvature scale but rather to the discrete quantum- 
gravity scale related to the Planck length, there is no a-priori reason why they should be 
small at low curvature. They can play a role in standard cosmological scenarios, for instance 
during inflation. Indeed, in such a combined scenario, the window allowed for inverse-triad 
effects is much smaller than the one for curvature or holonomy modifications. For inverse- 
triad corrections, a parameter range of about four orders of magnitude is consistent with 
observations [59l 1146] . while curvature corrections have an allowed range of about ten 
orders of magnitude, one compared to the ratio of densities in observationally accessible 
regimes to the Planck density. There are also indications of interesting and characteristic 
effects in non-Gaussianity [147] . although the required equations of motion second order 
in inhomogeneity still have to be made consistent. By inverse-triad effects, loop quantum 
gravity becomes falsifiable. 

For a more detailed review of phenomenological implications, see |148j . 

7 Outlook 

Any quantum system can be evaluated consistently and reliably only when all possible 
quantum effects and the relevant degrees of freedom are taken into account. For quantum 
cosmology, this means that one must go beyond the traditional minisuperspace models and 
find consistent extensions to inhomogeneity. Quantum-representation issues then become 
much more involved, but can be handled for instance with methods of loop quantum grav- 
ity. In this canonical setting, a large set of effective techniques, described in the main part 
of this review, is now available. These methods allow one to forgo ad-hoc assumptions, to 
implement full (but possibly deformed) space-time covariance, and to derive a complete 
phenomenological setting in which all relevant quantum effects are included. Unlike in tra- 
ditional canonical quantizations and derivations of wave functions, there do not appear to 
be major obstacles on the way toward systematic comparisons with observations. Quantum 
cosmology is therefore empirically testable. 

Concretely working out all terms and studying the necessary parameterizations of quan- 
tization and state ambiguities still remains to be completed. Even in isotropic models be- 
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yond the harmonic one, not much is known about the evolution of generic quantum states 
and the robustness of singularity avoidance (see Sections 12 .21 and 14")) . Control on inhomoge- 
neous modes, necessary for most physical questions in cosmology and an understanding of 
quantum space-time, remains poor in strong quantum regimes, suffering from quantization 
ambiguities and the difficult anomaly problem (see Sections 12.11 and 151. Effective tech- 
niques, especially effective constraints, have relieved some of the pressure caused by the 
failure of traditional methods to address these problems, but they have not been evaluated 
in sufficient detail to provide a reliable view on Planckian stages in cosmology. Further 
in-depth investigations are required to change this situation and to provide a complete and 
reliable phenomenology of quantum cosmology. 
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